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Abstract 

In this review article we discuss the recent progress in studying baUistic transport for charge 
carriers in graphene through highly inhomogenous magnetic field known as magnetic barrier in 
combination with gate voltage induced electrostatic potential. Starting with cases for a single or 
double magnetic barrier we also review the progress in understanding electron transport through 
the superlattices created out of such electromagnetic potential barriers and discuss the possibility 
of experimental realization of such systems. The emphasis is particularly on the analogy of such 
transport with propagation of light wave through medium with alternating dielectric constant. 
In that direction we discuss electron analogue of optical phenomena like fabry perot resonances, 
negative refraction, Goos-Hanchen effect, beam collimation in such systems and explain how such 
analogy is going to be useful for device generation. The resulting modification of band structure of 
dirac fermions, the emergence of additional dirac points was also discussed accompanied by brief 
section on the interconvertibility of electric and magnetic field for relativistic dirac fermions. We 
also discuss the effect of such electromagnetic potential barrier on bilayer graphene in a similar 
framework. 
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I. INTRODUCTION 



Physically Maxwell's equation that describes the propagation of electromagnetic 
wave/light through free space and dielectric medium, namely 

and Schrodinger equation that describes the time evolution of probability amplitude or de 
Broglie wave in quantum systems, namely 

correspond to the different type of phenomena. However, their striking mathematical sim- 
ilarity indicates that a large number of wave like phenomena will occur in either of these 
cases. Whereas the motion of a isolated single electron in presence of potential barrier is 
described by the wave Eq.([2]), in real materials due to the scattering by other electrons, 
impurities etc. the motion of the charge carriers is generally diffusive. Therefore, such wave 
propagation based description of the electron transport becomes only meaningful if the elec- 
tron mean free path is of the order the typical sample size. Such a transport regime is called 
the ballistic transport regime and in that regime the similarity between the propagations of 
transport electrons and electromagnetic wave or light promises rich dividend [III2]- 
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A particular consequence of the similarity between the wave equations ([T]) and ([2]) is 
that for a monochromatic wave with a given frequency in Eq. [T] the role played by the 
dielectric constant of a given medium has its corresponding analogue in the the potential 
landscape for a stationary solution of the Schrodinger electron described by Eq. [2j Thus 
some of the effects that one obtains by spatially modulating the potential in a Schrodinger 
equation can be reproduced for the light also by spatially modulating the dielectric constant 
of the medium. Indeed this was pointed out in seventies in the pioneering work by Yariv and 
collaborators [3] that the propagation of light through a medium with periodically modulated 
dielectric constant will lead to the similar band structure of transport electron as observed 
in Kronig-Penny model [1]. In this context it is particularly useful to point out that for 
conventional two dimensional electron gas (2DEG), the analogy between transmission of de 
Broglie waves satisfying the Schrodinger equation through a one- dimensional electrostatic 
potential O |6] and light propagation in linear dilelectric medium is well established and led 
to the development of a number of applications. 

The recent discovery of graphene [7HTU] added a new twist to this well established optical 
analogy of ballistic electron transport and vice versa. In graphene, electrons near the Fermi 
level, namely the transport electrons no more obeys quadratic dispersion law, a typical char- 
acteristics of their non-relativistic nature, but rather obeys a linear dispersion relation at 
or near the Fermi surface, an archetypical of the ultra relativistic massless particles. Hence 
they are called massless Dirac fermions. It has now been established that transport of such 
massless Dirac fermions in the presence of an electrostatic potential barrier is analogous to 
negative refraction through metamaterials [TlHT3] . The relativistic behaviour of graphene 
electrons also leads to Klein tunnelling [13], where a massless relativistic particle can tun- 
nel through any potential barrier above the Fermi surface, invalidating the possibility of 
confining it using such potential barriers. 

In the context of electronic transport in graphene devices that are proposed for graphene 
based electronic in future, the local carrier concentration is controlled by one or more local 
gates. Particularly, in the cleanest devices which satisfy the criterion of ballistic transport, 
transport signatures of the relativistic nature of the charge carriers of graphene are observed. 
A vast body of theoretical work ( for a review see [15]) has already been devoted to the study 
of multigated graphene superlattices which can form a basis for a new kind of electronic 
optics based on graphene as an electronic metamaterial. However the inability of confining 
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the electrons using a electrostatic potential barrier severely limits the applicability of such 
devices. A solution to that problem was proposed by De Martino et al. (TU] when it 
was pointed out that well localized magnetic field dubbed as magnetic barrier can confine 
massless dirac fermions in graphene. This suggests that one way of making high mobility 
graphene based electronic devices is to make locally gated structures in conjunction with 
such magnetic barriers. Given this context it is extremely important to see what will happen 
to the above mentioned optical analogy when it applies to the cases of such magnetic barriers. 

However it may be pointed out that the study of transport properties of two dimensional 
non relativistic electron system in the presence of a transverse magnetic field remains mostly 
confined to the case where the magnetic field is uniform on the scale of sample size through 
phenomena such as Integer and Fractional Quantum Hall effect. Study of electron transport 
in inhomogenous magnetic field is relatively lesser known field [T7], even though inhomoge- 
nous magnetic field is being used for many other purpose for a long time, such as one in 
the famous Stern Gerlach experiment [TH]. In a homogenous magnetic field, the electrons 
execute cyclotron motion, the direction of the wavevector continuously changes, hence ren- 
dering any analogy with the monochromatic light propagation in dielectric medium is very 
unlikely. In quantum mechanical language the electronic wavefunction is localized over the 
scale of the magnetic length. Thus a direct analogy with a propagating wave is not possible. 

This situation however changed, when it was pointed out [19] that an optical analogy 
can be constructed for electron transport in graphene in presence of highly inhomogenous 
magnetic field [20] where one does not get bound state solutions, but rather scattering state 
solutions. This suggests that inhomogenous magnetic field not only confine massless dirac 
fermions [16], but the ballistic transport through such barriers can also be understood in 
terms of suitable optical analogy. Ref. [19] also showed that the analogy with light propa- 
gation in medium with periodically modulated dielectric constant can again be constructed 
by considering the transport of massless dirac fermions in graphene through the periodic 
arrangement of such highly inhomogenous magnetic field. Subsequent work [21] showed that 
a very rich transport regime can be realized for such graphene electrons by considering struc- 
tures where a local gate voltage is used simultaneously with such magnetic barrier. When 
such configurations are generalized to a superlattice structure, that significantly modifies 
the band structure of graphene electrons. Subsequently a large volume of theoretical work 
pointed out the relevance of such optical analogy and proposed device structures based on 
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that. Other properties such as Goos Hanchen shift of graphene electrons by these type of 
barriers, coUimation of the electron beam, generation of additional Dirac points was also 
discussed in further works. It was also extended in the case of bilayer graphene. In the 
current review article we plan to review these developments. The purpose is two fold. One, 
is to review the general theoretical framework in a pedagogical way, highlighting a number 
of significant work done in this subfield of graphene research. The other purpose is to indi- 
cate in very clear terms the possibility of designing interesting experiment and novel device 
structure that holds rich promises for graphene based electronics. 

The organization of the article is as follows. After giving a brief review of low energy 
description of charge carriers in graphene, first we explain how electron transport through 
scalar electrostatic barriers and magnetic barriers is explained in the language of geometrical 
optics. Next, we explain how such transport gets modified when various types of potential 
barriers are created atop magnetic barriers. We discuss Goos-Hanchen shift and fabry perot 
resonances in this type of structure. Then we provide an analysis of transport through 
an infinite series of such barriers pointing towards an effective way of changing the band- 
structure. Here we also discuss the collimation of the electron beam, emergence of extra 
dirac point etc. We then briefly discuss the issue of interconvertibilty of electrostatic poten- 
tial and vector potential due to magnetic field citing relativistic invariance of the equation 
obeyed by transport electrons in graphene. Then we come to the issue of electron transport 
through graphene bilayer in presence of such potential barriers. The dispersion relation for 
graphene bilayer is very different from that of monolayer graphene. The charge carriers in 
bilayer graphene have a parabolic energy spectrum, which means they are massive, similar 
to conventional non relativistic electrons. On the other hand, due to the crystal structure 
of graphene consisting of two sublattices, these particles are described by spinor wave func- 
tions similar to that of monolayer graphene. We study how the transport of such electrons 
in graphene bilayer is affected, again when exposed to combination of magnetic barriers and 
voltages. 

A. Low energy description of charge carriers in graphene 

The interesting band structure of graphene was pointed out by Phil Wallace [7] way back 
in 1947. After its experimental discovery by Geim, Novoselov, Kim and their collaborators 
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there are many excellent review articles (for example see [22] ) where the band structure, 
the resulting dirac fermion nature of the charge carriers are detailed. Here we provide a 
brief introduction to the low energy description of these charge carriers for completeness. 
Graphene is a one atom thick planar sheet of sp^ bonded carbon atoms. The carbon atoms 
are distributed at the edges of regular hexagons forming a honeycomb lattice. However, the 
honeycomb lattice is not a Bravais lattice and, from a crystallographic point of view, has to 
be described by two inter penetrating triangular lattices with two atoms (sublattice A and 
B) per unit cell. 

Each carbon atom has four valence electrons, three of which hybridize to form sp^ hy- 
bridized orbitals to form sigma bonds with neighbouring atoms. However, the fourth un- 
bound electron lies in the 2pz orbital which extends vertically above and below the plane. 
It is this electron which interacts with the periodic field of the hexagonal crystal lattice of 
graphene. If we denote the creation operator by a\{h\) for an atom on the A(B) sublattice, 
then the nearest neighbour tightbinding Hamiltonian has the simple form [231 [21] 

n = -ti ^{albj + H.c.) (3) 

where ti ^ 2.8meV is the nearest neighbour hopping parameter. The tightbinding eigen 
functions has the form of bispinor, whose components corresponds to the amplitudes on 
sublattice A and B respectively within the unit cell. This also means that in addition to its 
usual spin 1/2, electrons carriers an additional pseudo-spin 1/2 associated with its sublattice 
degree of freedom. The energy levels of the lattice can be determined by operating with the 
Hamiltonian on its eigen function to give the corresponding band structure: 

E = ±Uy/3 + f{k) (4) 

where f{k) = 2 cos(gj,av^) +4cos(gx3a/2) cos(gj,av^/2) with lattice parameter, a ~ 2.46^4°. 

In Eq.Q the plus sign refers to the upper (tt*) band and the minus sign refers to the 
lower (vr) band. The unique feature of the band structure of graphene [23 [2S] is that its 
valence band (vr) and conduction band (vr*) touch each other at 6 specific points, which are 
also the vertices of the hexagonal Brillouin zone of this honeycomb lattice ( see Fig. Out 
of these six points, only two K and K' are non-equivalent and all the rest are related to 
these two points by symmetry. These points are also known as the K points in the literature 
and their positions in the reciprocal space are: 
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(a) (b) 




FIG. 1: (a) Honey comb lattice of graphene showing two inequivalent carbon atoms (with different 
colour) in a unit cell (b) Band structure in honey comb lattice, with zoom in of the energy bands 
close to one of the Dirac points. 



Under ambient conditions, also known in the literature as zero biased doping, the Fermi 
level coincide with these corners of the hexagonal Brillouin zone. Thus, to understand the 
low energy charge transport properties of graphene one expands the full band structure close 
to the K point by writing q = K + k, with k <^ K and then by keeping the first order 
term and neglecting all higher order terms in A;, the resulting dispersion relation is of the 
form E = ±hvp\k\ such that low energy quasiparticles of graphene are described by the 
Dirac-like Hamiltonian [27H29] 



where vp ~ 10^ms~^ is the fermi velocity, cr = {a^jCy) is the vector comprising of two 
component Pauli matrices, and p = {px,Py) is momentum vector in x — y plane. Thus, 
charge carriers in graphene behave like massless relativistic fermions dubbed as massless 
Dirac fermions with only the velocity of light is replaced by the Fermi velocity. This way, 
the behaviour of transport electrons in graphene is very different from those in ordinary 
semiconductors where they have a parabolic dispersion like a non-relativistic free particle. 





(5) 
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In the following sections we shall describe how such charge carriers with ultra relativistic 
dispersion laws transports in presence of highly inhomogenous magnetic field dubbed as 
magnetic barriers in combination with local gate voltages. 

II. ELECTRON TRANSPORT IN GRAPHENE : OPTICAL ANALOGY 

A. Electron optics in the presence of scalar potential barriers 

In a two dimensional electron gas (2DEG) it is well-established that transmission of de 
Broglie waves satisfying the Schrddinger equation through a one- dimensional electrostatic 
potential is similar to light propagation through a refractive medium. Such transport can be 
understood in terms of phenomena like reflection, refraction and transmission, leading to an 
analogy between electron transport and light propagation [H El E] • When a non-relativistic 
electron in a 2DEG with quadratic dispersion, at Fermi energy E, is incident on a potential 
barrier V, its momentum parallel to the interface outside and inside the barrier is conserved; 
i.e., |pi|sin6'i = |p2|sin^2, where pi^2 are the momenta and 612 are the angles they make 
with surface normals respectively in the regions without and with potential barrier. This 
leads to the following Snell's law: 



Graphene charge carriers however do not have the quadratic dispersion, but instead 
behave as massless Dirac-Weyl fermions leading to a different set of transport phenomena 
P 4T01 [22l [30l [3T] . The optical analogues of such electron transport can again be constructed 
by considering the charge carriers in MLG incident on an electrostatic potential barrier 
V. Such charge carriers in monolayer graphene (MLG) obey the following Dirac-Weyl like 
equation: 



where = Vq in the barrier region < d) and vanishes outside, and = [\E'i, \I'2]"^ 

is the two component wavefunction with T denoting the transpose of the row vector. We 
assume V = V{x) to be the one dimensional potential. Due to the y invariance of potential 
barrier we consider the solutions of the form = il){x)e^^y'^ . On substituting this in 




(6) 




(7) 



9 



Eq.Q, two coupled equations in ipi and ip2 are obtained as: 

-I ^±ky ^2,1 = {E- Vi)ilj 



1,2 



The above two off-diagonal equations can be decoupled in terms of 'ipi^2 to give a 
Sclirodinger like equation of the form 



'2 



^1,2 = {E- Vi)\i^2 (8) 



which admits exponential solutions, namely ipi^2 oc exp(— i/c^.x). 

By parameterizing the energy momentum relation k'^.+ky = {E—ViY polar coordinates, 
we obtain 

ky = i?sin0, fcj.^ = Ecos(f), \x\ > d, 
ky = {E- V) sin 6, k^^ = {E -V) cos 6, \x\ < d. 

The momentum component along a straight interface should be conserved. Accordingly, 
using the electron momentum conservation in the y direction at the left interface x = —d, 
we obtain the Snell's law in the following form: 

{E -V)sme = E sin 

Clearly for E < V, the barrier acts like a medium with negative refractive index [11]. This 
is shown in Fig|2] 

sin 9 E 

n = = r (9) 

sin0 \E-V\ ^ ^ 

This corresponds to the electronic analogue of the well known phenomenon of negative re- 
fraction, which occurs in left handed metamaterials and which was first proposed by Veselago 
[12] and subsequently developed in more detail by Pendri and collaborators [13j. If one uses 
so called split gate voltage by placing the region V < E and V > E side by side ( Fig. |2]), 
the Fermi level can be tuned below and above the Dirac point. Since the region below and 
above the charge neutral Dirac points are respectively hole and electron states for such Dirac 
fermions, this creates p and n type region in graphene simply by changing the height of of 
the electrostatic potential barrier in a given region [TT]. The resulting structure is called 
graphene p — n junction. Consequently, a graphene p — n junction can be used as a Veselago 
lense for electron focussing. Similarly circular graphene p-n junctions has also been studied 
showing the formation of caustics [32] . 
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FIG. 2: Schematic of the transmission through a potential step for an electron with (a) E > V (b) 
E <V. 

B. Klein Tunneling 

While electrostatic potential barriers can indeed manipulate charge transport in graphene, 
an electron could tunnel through any high barrier in contrast to the conventional tunnelling 
of non-relativistic electrons [HI [33]. This behaviour, called Klein tunnelling in graphene, 
leads to several observable transport effects related to transport some of which have been 
demonstrated in graphene [311 ES] and also in carbon nanotubes [36]. Below we review 
briefly the work on such tunnelling problem for graphene charge carriers in the presence of 
an electrostatic potential barrier. This analysis will then be extended to similar transport 
through inhomogenous magnetic fields and will be followed by comparing and contrasting 
these transport to electrostatic and magnetic barriers, which is the central topic of this 
review article. 

The wavefunction solutions in the presence of a potential barrier can be obtained by 
solving EqjSj These solutions in any region of space can be written in terms of linear 
combination of forward and backward propagating plane waves such that 



Q^^iqxx _|_ fj^-iq^x^ \x\ < d 

te*''^^, X > d 
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^2(3;) = {s' (ae^«-^+*^ + 6e-^^-^-*^) , \x\ < d (10) 

where fci? = 27r/A is the Fermi wavevector, = kpcoscj) and ky = kpsiiKp are the 
wavevector components outside the barrier, = ^ {E — Vy/h?Vp — k^ ,9 = tan~^{ky/qx) 
is the refraction angle, s = sgn{E) and s' = sgn{E — Vq). 

By taking into account the continuity of wave-function components ipi ^ind 1^2 at the 

boundaries of the barrier the transmission and reflection coefficients are obtained as 

2ss'e~*''^^ cos0cos6' 

~ ss'[e-'i-^ cos((/) + 6') + e'l-^ cos{(j) - 9)] - 2i sin q^D ^ ' 

2ie*'^ sin ga;D[sin — ss' sin 6*] , , 

— (12j 



ss'[e-^i-^ cos(0 + 9) + e'l-^ cos{(j) - 9)] - 2i sin q^D 
Here D = 2d is the width of the barrier. Clearly q^D = mr, n = 0, ±1, ... corresponds 
to the resonance condition at which the barrier transmission is unity. This corresponds to 
usual resonant tunnelling through a potential barrier and occurs for non relativistic electron 
as well. However remarkably, the above expression also shows that barrier always remains 
transparent for normally incident carriers, ie. T = 1 when = 0. This feature of anomalous 
tunnelling in graphene is peculiar to the Dirac like spectrum for graphene charge carriers 
and is explained as follows: As explained briefly in section hK the gapless, conical spectrum 



of graphene is the result of intersection of two cosine like energy bands originating form 
sublattices A and B. (represented in green and red for A and B sublattice respectively, see 
figjl|. Due to this, the sublattice degree of freedom i.e. pseudospin should remain fixed on 
each branch. This also means that an electron with energy E and possessing wavevector k 
will originate from the same branch as a hole with energy —E and possessing wavevector 
— k. Above two features allows the introduction of chirality for graphene charge carriers. 
The term (cr ■ p) present in the Hamiltonian ^ gives the chirality, namely projection of 
pseudospin on the direction of motion, 

\P\ 

and has two possible eigen values, +1 (for electrons) and —1 (for holes). This operator is 
same as the usual helicity operator for 3 + 1-dimensional relativistic electrons that obey Dirac 
equation [37], but redefined for charge carriers in graphene which obeys 2 + 1 dimensional 
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Dirac-Weyl equation. In the absence of any potential, the chirahty operator commutes with 
the Hamiltonian and is therefore a conserved quantity. For normally incident carriers, the 
chirality remains conserved even in the presence of an external electrostatic potential V{x)l. 
This can be shown as follows: For normally incident carriers, 

Chirality, C ■ 



Px 



Hence 



Also the velocity operator 



[C,H] = [a,,k,a, + V(x)l]=0 



Vx = -i[x,H] 



Thus the velocity operator is same as the chirality which is a conserved quantity for 
this case [211 |38]. The velocity along the x direction is therefore a constant of motion 
and thus cannot be reversed. This leads to perfect transmission through such a barrier 
at normal incidence. This is what describes Klein tunnelling in MLG [E]. This is to 
be noted here that a rectangular electrostatic barrier in graphene assumes that the gate 
voltage-induced doping changes abruptly at the edges . In reality, however, the doping level 
continuously changes and thus the edge of the potential barrier is actually smooth and not 
sharp. This issue has been considered for scalar potential barrier in Ref. [39], where it was 
found that a potential which is smooth on the scale of the Fermi wave length for small angles 
of incidence, T(0) = exp^—irkpD siii^ (j)) where D is the barrier width. A comparison with 
the transmission expression given in EqjlT] shows that the assumption of the rectangular 
barrier captures the effect of Klein tunnelling correctly. However, at other angles close to 
the normal incidence it overestimates the transmission. Transmission through trapezoidal 
barrier was also analysed in Ref. [IQ] which combines the effect of a smooth barrier and a 
rectangular barrier. 



On expanding (j) in eq.(12), close to the normal incidence, namely = 0, by substituting 
(f) = 6(j)o, it can be seen that the reflection amplitude r undergoes a tt phase jump when the 
incident angle goes from positive to negative value . At zero magnetic field, at non-normal 
incidences, the two consecutive reflections on the two p-n interfaces occur with opposite 



angles 6i = 6 and 62 = —9 (Fig, 20). When such a system is placed in a transverse magnetic 
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field B = Bz, the electronic trajectories bend in the presence of magnetic field. And above a 
critical field value B > Be, trajectory bending becomes sufficient to make the two consecutive 
reflections occur with the same incident angle 6i = 62 which is exactly the case what happens 
at normal incidence in the absence of magnetic field. This suddenly adds vr to the phase 
accumulated by an electron between two reflections and shifts the interference fringes by 
half a period [H]. The observation of this half-period shift in Fabry Perot interference 
fringes is therefore a direct evidence of perfect tunnelling at normal incidence |3l]. We 



analyse this in more detail in section IV D An alternative way of demonstration of Klein 
tunnelling in graphene p-n junctions has been addressed in [33]. This was done by probing 
the transition from clean to disordered transport across a single steep p-n junction. Very 
recently, an angle dependent carrier transmission probability in graphene p-n junctions has 
also been investigated experimentally [12] and theoretically |13] where it is shown that chiral 
tunnelling can be directly observed from the junction resistance of a tilted interface probed 
with separate split gates. 



III. ELECTRON TRANSPORT IN THE PRESENCE OF INHOMOGENOUS 
MAGNETIC FIELD PROFILE 

While the above mentioned absence of backscattering by a scalar potential due to Klein 
tunnelling is a very interesting phenomenon, this implies that confining transport electrons in 
graphene by a potential barrier is not possible in a conventional way. This leads to problem 
in device making. For example, even though it is easy to make a p — n junction in graphene, 
reverse biasing such junction will be very difficult for Klein tunnelling. For the designing of 
graphene based electronics it is crucial to attain confinement of electrons within a mesoscopic 
or nanoscopic size of the sample. For this reason, several alternatives have been suggested. 
One way is to exploit the fact that suitable transverse states in a graphene strip may allow 
one to circumvent Klein tunnelling [H]. Other schemes that have been proposed include 
gated nanoribbons [45J, gated |16J or doped ^47j bilayer etc. Yet another possibility was 
demonstrated theoretically in [IB], [IE] by making use of external inhomogeneous magnetic 
fields applied perpendicular to the graphene plane. This brings us to the the main focus 
of this review article- to investigate the electron transport in the presence of inhomogenous 
magnetic field viz. magnetic barrier(s) in graphene. Before proceeding further we first 
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FIG. 3: Magnetic field under the stripe corresponding to the (a) perpendicular (b) parallel config- 
urations. The magnetic field is given at the following distances from the magnetic stripe: zq = 0.15 
(solid curve), zq = 0.3 (dashed curve), and zq = 0.5 (dotted curve). 

discuss the experimental strategies for creating such inhomogenous magnetic fields and point 
out that such magnetic barriers are already in very much use in other materials. 

A. Producing inhomogenous magnetic field profile: Experimental strategies 

A microscopically inhomogenous magnetic field can be created for charge carriers in 
graphene by placing a graphene sheet in close proximity to long magnetic stripes that 
produce highly localized magnetic fields. Such field profiles can be generated using de- 
magnetizing fields produced at the edges of narrow stripes made with hard ferromagnetic 
(FM) materials of either perpendicular or in-plane anisotropy. The magnetic field of such 
structures is given as 

B = B{x, Zq)z = Bq[K{x + d, Zq) — K{x — d, Zq)]z 

where K{x,zo) = — Jf^^i for perpendicular magnetization, and K{x^zq) = — J^Jj^a for mag- 
netization parallel to the graphene sheet. Bq is a constant dependent on the aspect ratio of 
the stripe. 

For a given value of Zq we plot in Fig|3|^a), (b) the profile of such a magnetic field. The 
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major component of the demagnetizing field reaches the graphene monolayer and produces 
the desired field profile. Though there is always some component of the demagnetizing field 
that will give rise to undesired fringe fields in other directions, these can be substantially 
lowered by suitable magnetic design of the stripes. Such nanostructures are routinely used 
in magnetic recording media [19]. Materials such as CoCrPt produce fields of 1 Tesla close 
to the surface with bit lengths ranging from 50-100 nm. For example, in [50], isolated tracks 
of single-domain magnetic islands have been fabricated using focused ion-beam lithography 
for both perpendicular and parallel anisotropy using CoCrPt. Here, patterns with successive 
magnetizations pointed along opposite directions were achieved. In another recent work, off- 
axis electron holography has been used to probe the magnetization structure in high density 
recording medium by using perpendicular magnetic anisotropic (PMA) recording medium 
[51] . The direct imaging of magnetization done shows that the foils of PMA material consist 
of successively reversed highly stable domain structures of few ten's of nanometer size. In 
practice, one can also change the strength of the magnetic field by suitable adjusting the 
width of such PMA material. Precise design of read-write structures for recording individual 
bits at these dimensions has also been achieved [52]. Using the above discussed techniques, 
typical magnetic barriers can be patterned down to 50-100 nm widths. 

Apart from the above mentioned method of using ferromagnetic materials [53] there are 
several other ways in which inhomogenous magnetic fields have been experimentally realised 
in conventional semiconductors. One such way is through the integration of superconducting 
elements, which may be used to screen the externally applied magnetic field in accordance 
with Meissner effect. When an external magnetic field is applied the flux lines will be ex- 
pelled from the superconductor due to the Meissner effect. If the latter is close enough 
this will result in an inhomogeneous magnetic field in the 2DEG [5l][55] . This is shown 
schematically in [Figj3](d)]. As an example, in ref. [54j, it was shown that with type II super- 
conducting films deposited on top of the two-dimensional electron gas in a GaAs/AlGaAs 
heterostructure, the distribution at the 2DEG takes the form of flux tubes which are much 
narrower than an electron inelastic scattering length. The superconducting materials used 
in these experiments were typically 200nm of lead (Pb), 400 nm of a lead/indium [Pb(l at. 
% In)] alloy, and 200 nm of niobium nitride (NbN), as being materials whose characteristic 
superconducting length scales span a broad range. Growing the heterojunction on top of a 
pre-etched (nonplanar) substrate can also give rise to inhomogeneous magnetic fields [5S] - 
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FIG. 4: (color online){a) Cobalt finger gates at the surface of a 2DEG. (b) Magnetic modulation 
profiles at different tilt angles of the magnetization, (c) MFM image of the cobalt grating magnetized 
along the short axis of the stripes, in the plane of the 2DEG. The magnetic poles appear as the 
black lines, (d) MFM image of the grating magnetized along the long axis. The magnetic poles 
are absent. Parameters: a = 4QQnm, d = 200nm, h = 160nm and zq = 90nm. [This figure is taken 
from [17j, printed with permission.] 

|58j . When a uniform external field is applied on such a system, the angle between the field 
direction and the normal to the 2DEG depends on the tilt of the facet and therefore will 
have a normal component of magnetic field which varies spatially across the sample. Since 
it is only the normal component of magnetic field that influences the transport, by varying 
the facet length and angle, as well as the angle between the applied field and the normal to 
the substrate, a wide variety of field profiles can be generated. 

In other studies |59]-|6l] a similar yet simpler approach was used based on the idea that 
if a thin magnetic film is placed on top of a heterostructure its in-plane magnetization can 
be saturated in an external in-plane magnetic field [Figjsj^c)]. In this case the out-of-plane 
component of the fringe field under the edge of the film creates a magnetic barrier for electron 
transport with a width of the order of 100 nm determined by the separation of the 2DEG 
from the magnetic film. Magnetic field strengths of more than 0.5 T have been realized in 
this way A detailed discussion on experimental methods used for applying magnetic 



17 



5 40 80 5 40 80 

E(meV) E(meV) 



FIG. 5: Transmission in the presence of finite width magnetic barrier in monolayer graphene, 
D = 2d = 1.51b {si)B = 0.5T (b) B = 3T. 

modulations to 2DEGs can be found in the review article |17j . 

B. Magnetic confinement of massless dirac fermions in graphene 

In the previous section we reviewed the experimental progress of realizing inhomogenous 
magnetic field or magnetic barrier over a wide range of length scales in a number of sys- 
tems. Mesoscopic transport in presence of inhomogenous magnetic fields has been studied 
theoretically [62] , [63] - [6^ for non relativistic electrons . Related experimental studies also 
took place in conventional semiconductor heterostructures, e.g., transport in the presence of 
magnetic barriers [65] and superlattices [66], magnetic edge states close to a magnetic step 
[67j . and magnetically confined quantum dots or antidots [68]. In the following section B we 
review the transport of massless Dirac fermions in graphene in the presence of of magnetic 
barriers with some suitable examples. 

The magnetic barrier profile chosen here is in the form of step like function [T6l HE] as 



, and the vector potential chosen in Landau gauge Ay{x) = [—Bd, Bx, Bd] in region x < —d, 
\x\ < d and x > d respectively. 

Typically, the equation of motion of charge carriers at the two inequivalent Dirac points 
K and K' are decoupled from each other [69] in the absence of intervalley scattering which 
gives a wave vector change of the order of 2kF. We assume for the cases to be considered 




(13) 
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this condition is satisfied. Under this situation motion near each Dirac point can be treated 
separately. The stationary solutions near each such point is given by a two component spinor 
and satisfies the equation 

Vpi^a^ T ^7ry)^2,l = ^^1,2 (14) 

Here, vp is the fermi velocity {vp ~ c/300) and tt = p + ^ A. Using Ib as the unit of length 
scale such that x = [|-and ^ as the unit energy such that e = and ^{x,y) = il)[x)e^^yy 
in the Landau gauge, Eq|l4]can be rewritten as 



d 



± kJ 



which can be decoupled to obtain the following Schrddinger like equation: 



Ay{x) 



^2.1 = e^i,: 



(15) 



Q2 

+ {kylB+xf + 1 



dx 

d 



^1,2 = eVi,2 < d 



dx^ 



— + {kyls =F d/lsY 



^1,2 = e 



1,2 



(16) 



where (— /+) sign in the above expression corresponds to region {x < —d)/{x > d) respec- 
tively. 

Upon solving, we obtain the solutions outside the barrier as propagating solutions while 
inside the barrier the solutions are localised and can be written in terms of parabolic cylin- 
drical functions [70] . 



^i(x) 



X < —d 
aD ^2 /2^i{y/2{x + kylB)) 
+bD ^2 /2^i{-V2{x + kylB))-, \x\ < d 



te^^'^^ , X > d 



(17) 



(18) 



a^D,2/2(v^(x + kyh)) 
Vh^D,2/2{-y/2{x + kylB)), \x\ < d 

where kr = kp COS (p^ — COS0', D stands for parabolic cylindrical function, and 
and 0' corresponds to the incidence and emergence angle in regioin I and III respectively. 
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The emergence angle can be obtained by conservation condition of ky as: 

2d . . 

sin0 = — h sin0 (19) 

An alternative way of writing the above equation is 

$ 

e sin 6' = 21: he sin 6 

$0 

where $o = hc/e and $ = 2BdlB, is the net flux enclosed by an area 2dlB- The second 
form particularly shows how the flux enclosed by the area containing the barrier changes 
the refractive index of this region in a non specular way. 

The transmission through such magnetic barrier can be obtained by matching the wave- 
function at the boundaries of the barrier, T = |f|^P, where the prefactor ensures the current 
conservation in the system. The transmission as function of incidence angle and incident 
energy E is shown in figjSj As one can see, Eq. 17 and [18] depicts the formation of lo- 
calised states inside the barrier region. However these states are bound only in the direction 
perpendicular to the barrier but can propagate in the direction along the inhomogeneity of 
the magnetic field. The transmission, for this reason, depends very strongly on the incident 
wavevector and facilitates the possibility to construct wave- vector filter Fig. |5] depicts 
this situation clearly. The requirement of the propagating solutions in the incident region 
and region of exit determines the range — 7r/2 < < sin^^ ^1 — ^ j beyond which no trans- 
mission is possible, thereby leading to the possibility of confinement. This is how magnetic 
barrier turn reflective for a set of wave vector for the massless Dirac fermions. In the follow- 
ing section we shall illustrate this property more clearly by building clear optical analogy 
of electron transmission through the extremely inhomogeneous magnetic barrier where the 
magnetic field profile can be approximated as a delta function. 



IV. MASSLESS DIRAC FERMIONS IN INHOMOGENOUS MAGNETIC FIELDS: 
AN ANALOGY WITH LIGHT PROPAGATION 

A. Magnetic vector potential (MVP) barrier 

The issue addressed here is whether the transport through such magnetic barrier can be 
understood in terms of propagation of light through medium with changing dilelectric con- 
stant using well known ideas in geometrical optics[19]. Even though such idea was developed 
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for the ballistic electron transport through electrostatic potential barrier, generalizing this 
concept in presence of a magnetic field was non trivial since in the presence of magnetic field 
electron executes cylotron motion and the wave vector k is not a good quantum number for 
such motion. Thus to yield a Snell's law like one in Eq. ^ in presence of a magnetic barrier, 
it is essential that the electron should not complete its half cyclotron radius beyond which it 
get completely localised inside the barrier. This situation can be realized in the presence of 
a highly inhomogeneous magnetic field. Particularly, if the range of inhomogeneity is much 
smaller than the cyclotron radius, one is left with plane wave-like scattering states. Thus, 
such field profile scatter the electrons in the same way as an electrostatic potential does. 

For this to be valid, two conditions must be satisfied. First, the magnetic length Ib = 
\J~^ should be similar in order to the width of such magnetic barriers. Second, the de 
Broglie wavelength should be much larger than the barrier width so that the electron 
will not see the variation in the vector potential inside the barrier. An extreme case of 
such an inhomogeneous magnetic field is singular magnetic barrier - the one introduced in 
[TSH2I1 EH mHT5| . which gives rise to step-function like magnetic vector potential (MVP) 
barriers. Here, we shall consider the scattering of massless Dirac fermions by such singular 



magnetic barriers as in Eq,20 Such field profiles having highly localized field variations are 



well known in literature and have been realized experimentally as we have already pointed 



out in the section III A Patterned stripes down to length scales as low as 10 nm have 
been realized using nanolithography [71]. It is possible to achieve field profiles at even 
smaller dimensions using domain walls with widths in the range of 10-50nm and magnetic 
nanostructures down to 5 nm [75l|76] and 0.15 nm [77] having highly localized field variations 
which can be approximated as such delta function barriers. 

1. Electron optics with MVP harrier 

The most elementary of such magnetic field profile that satisfy the flux line closure 
condition and its vector potential in the Landau gauge can be given as : 

B = BIb[5{x + d)- S{x - d)]z (20) 
Ay{x) = BlBe{d^-x^)y 
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FIG. 6: Monolayer graphene with ferromagnetic stripes having magnetizations perpendicular (a) 
and parallel (b) to plane. The magnetic field B (red) produces a magnetic vector potential A 
(blue). Single MVP barriers are formed in (a) and (b). 

In the absence of intervalley scattering, the charge carriers in graphene in the presence of 
above potential form can be described by EqjM] which can be rewritten in the form of the 
following coupled equation 

V^2,i = eV^i,2 (21) 

Here A = 1 for |x| < d and A = for > d. The coupled equations Eq|2T] can be 
decoupled easily and results in a Schrodinger hke equation of the form 

V'1,2 = (22) 

Thus, in the barrier region —d < x < d electrons see a momentum-dependent barrier of 
height [ky — sgn{e)-^]'^. Hence such a barrier will be termed as magnetic vector potential 
barrier (MVP). 

Since magnetic field doesn't do any work, the energy conservation gives k'^ + = k\ for 
|x| > c? and + {ky + = k\ ioi \x\ < d. By parameterizing the above two equations in 
polar co-ordinates as shown in Fig. [7], we obtain the relation 

sin \e\ = sin |0| + sgn{(f))-^, -n/2 < < 7r/2 (23) 

The situation is depicted in FigjT} For a wave incident with positive (j) wave vector bends 
away from the normal whereas for a wave incident with negative incidence angle the cor- 
responding wave vector bends towards the surface normal inside the barrier region. Here, 
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FIG. 7: Asymmetric refraction through a single MVP barrier. 



9 is the angle of refraction (denoted as 61^2 for positive and negative angles of incidence). 
Clearly, the Snell's law for electron waves in such magnetic barriers is not specular as it is 
for light waves on smooth surfaces or for the incidence of the electrons on an electrostatic 
potential barrier [TT], [H]. 

According to Eq|23| when \sin \9\\> 1^ 6 becomes imaginary and the wave in the second 
medium becomes evanescent. In the language of optics this corresponds to total internal 
reflection (TIR). According to Fig. [sj this will happen when sin|^| > 1 for < < | 
and when sin |6'|<— Ifor— 1<</)<0. In the latter case, this requires the wave vector 
to be negatively refracted |39j at sufficiently high magnetic field before TIR occurs. It also 
follows that for a given strength B the magnitude of critical angle of incidence \(j)\ = (pc for 
TIR is higher for — | < < as compared to the one for < < |. Because of TIR 
the transmission on both sides of Fig |8] drops to beyond a certain value of and for fixed 
B field this value is higher for negative angles of incidence . Thus we can understand the 



refiective nature of such barrier using the corresponding Snell's law given in 23 To illustrate 
this we shall now calculate the transmission through such barrier. 

2. MVP barrier: Transmission probability 



On solving Eq, 22 , the corresponding wavefunctions in any region of space can be written 



in terms of a linear superposition of forward and backward moving plane waves. Then 
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FIG. 8: Transmission T as a function of (p and E for single MVP barrier (a) B = 0.5T , d = 51b 
(b) B = 3T , d = 20lB. 

continuity of the wavefunction at the boundaries of the MVP barrier can be used to calculate 

the transmission coefficient as 

^_ 2ss'e~*''^'^ cos (/) cos 6* 

~ ss'[e-^i-^cos{4> + 9) + e*9-^ cos(0 - 6*)] -2i sin q^D ^ ' 

Here D = 2d. And s and s' are given by sgn(e) and are both +1 for electrons when 

only magnetic fields are present. It may be pointed out that the form of the transmission 



expression is same as the corresponding one for electrostatic potential barrier (11). But 
because the wave vector inside the barrier is different as compared to those in presence 
of scalar barrier, the transmission profile is fundamentally different. FigJs] shows how the 



Transmittance, T = t*t changes with the angle of incidence. As suggested by Eq,23, the 
transmission is clearly asymmetric unlike the case of electrostatic potential fU. 



For high electrostatic barriers such that V » Ep, the wavevector is given by 



^^KiJ^ which is real. The corresponding transmission for electrostatic potentials is 

y COS^ (t) 



1 — cos{qxD) sin^ 

and is 1 at = 0, exhibiting Klein tunnelling for massless Dirac fermions [H]. In contrast 

1 

kplB 



at high values of the magnetic barrier since t— i— oc vS, the wavevector given by 



k% — {ky + = —K^ < 0. This leads to TIR and not Klein tunnelling. As seen in Fig. 
the magnitude of critical angle beyond which TIR occurs is lower for a higher magnetic field. 
Thus, a stronger MVP barrier leads to higher refiections as opposed to complete transmission 
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at normal incidence by high electrostatic potential barrier. Complete transmission occurs 



only for q^D = nrr as given by Eq 24 This corresponds to resonant tunnelling for Dirac 
electrons and happens in the same way as for nonrelativistic electrons, appearing as a number 
of peaks in the plots in Fig|8] . The number of tunnelling peaks increases with barrier width 
for both MVP and electrostatic barriers. 

In conclusion, we see that the reflectance can be controlled by suitably modifying the 
strength and locations of the magnetic barriers and thereby changing the refractive index of 
the intervening medium in a novel manner. This principle could be the basis of structures 
like magnetic waveguide, where the reflection must be high at desired propagation angles, 
here, this could be manipulated by changing the magnetic field. Also for strutures like 
the resonant cavity, high reflection is needed near normal incidence. Geometries such as 
three-mirror or four-mirror cavities fi9\ could be used for high reflection at other angles . 



B. Combined MVP barrier and electrostatic potential (EMVP) barrier 

The charge carriers in graphene have a linear band structure albeit only close to the 
Dirac point, coincident with the Fermi level Ep in undoped graphene. Small electrostatic 
potentials greatly affect electron states by shifting the Dirac point with respect to Ep and 
causing the graphene sheet to behave as either an electron-deficit (p-type) or a hole-deficit 
(n-type) material. This situation was also experimentally verified by observing electron hole 
puddles [151 [78] near the Dirac point in monolayer graphene. Thus, the effect of electrostatic 
potentials on any proposed graphene structure must be included. 

Quantum hall effect in gate controlled p-n junction in graphene in presence of uniform 
magnetic field was already studied experimentally [7D] which reveals new quantum hall 
plateaus. The top gate geometry was utilized in controlling the edge channels in the quan- 
tum Hall regime and with control over local and global carrier density. The effect of local 
modulation of charge density and carrier type in graphene field effect transistors using a 
double top gate geometry has also been studied [80j. Also, in a recent experiment by S. 
Pisana et al.[HT], the enhanced magneto-resistance of a monolayer graphene sheet has been 
measured by connecting it to two voltage and two current terminals and simultaneously 
exposing it to various magnetic field strength at room temperature. The differential voltage 
as a function of the magnetic field has been plotted thereby analysing the joint effect of 
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magnetic field and the applied voltage on the magneto-transport particularly close to the 
Dirac point. From the analysis of the magneto resistance data it was inferred that the band 
structure of transport gets strongly modified in presence of voltage and magnetic field. The 
above experiment result clearly shows that simultaneous application of voltage and magnetic 
field strongly influences band structure, as we'll also see in the following section. For most 
of the above experiments the magnetic field applied is homogenous over the typical size of 
the sample. A complimentary case where as for our case we consider magnetic barrier which 
is a highly inhomogenous magnetic field over the typical size of the sample. 
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FIG. 9: EMVP barrier structures for graphene. In (a)-(c), the magnetic field is applied using 
patterned ferromagnetic (FM) lines with perpendicular anisotropy and two potentials +Vg and 
—Vg are applied by separate conductor lines. As given in the text, effective potential induced in 
the graphene is V. In (b), the magnetic strength of the middle line is doubled by either using a 
different FM material or by larger dimensions. In (d), the magnetic field is produced from the two 
edges of a FM stripe with in-plane anisotropy and the same stripe is also used to apply one of the 
potentials. 

We first discuss briefly some of the practical issues when a electrostatic gate potential is 
applied to the graphene sheet. A gate voltage ±Vg can be applied using a metal electrode 
and separating the electrode from the graphene layer with an insulating oxide layer. This 
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oxide layer which serves as a dielectric medium between graphene sheet and metal electrode 
can be made as thin as several nanometers. When such a voltage Vg is applied locally, it 
induces electron (hole) doping ±(7„ proportionally and changes the carrier concentration 
within the channel. This shifts the undoped Fermi level Ep from the Dirac point by an 
amount V = sgn{an)^/onhvF, where sgn{an) is the sign of the induced charge. Due to 
this a local potential barrier V is created which is equal to the difference between the local 
Fermi level and the Fermi level in the undoped region. Since the typical breakdown strength 
of dielectrics such as alumina and fused silica are around 10 — 20MV^/m, a graphene layer 
could be subjected to IV applied across a lOOnm thick dielectric. Following Refs.[SllS], gate 
voltages Vg of upto ±100\^ have already been applied to graphene flakes. The electrostatic 
gate potentials would be generally applied by separate conductors placed suitably in different 
planes than the FM stripes, but in some cases, the FM stripes could also be used to apply 
voltages. 

1. Electron optics with EM VP barrier 

Below we discuss the transport in the presence of a single MVP barrier and a commen- 
surate electrostatic step potential (EMVP) barrier. 



V{x) = V, Ay{x) = BIb, \x\ < d 

V{x) = 0,Ay{x) = 0,\x\>d (25) 

We assume ^ > and set the incident energy at E = Ep, where Ep is the Fermi level located 
at the charge-neutral Dirac point in monolayer graphene and the voltage V is measured with 
respect to Ep. The entire treatment that follows is applicable in the neighbourhood of the 
charge-neutral Dirac points as long as the dispersion remains linear. In addition, we also 
assume that scattering by an EMVP barrier is not strong enough to break the degeneracy 
of the K and K' points. In the region —d < x < d a.t either of the K points, the motion will 
be described by 

vJ ''-'A (A=Eh] (26) 

\7r^ + l7Ty V/Vp j \i)2) \V^2/ 
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FIG. 10: (a) Function Sp as a function of voltage (b) Ray diagram of electron propagation through 
a single EMVP barrier. With increasing V, the refracted angle 62 increases continuously in the 
clockwise direction. Dashed rays are negatively refracted waves. 

Here tz = p + ^A. In the Landau gauge, the stationary solutions can again be written as 



Substituting this in Eq,26, we get the following coupled equations 
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The above equation can be decoupled to yield 
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The corresponding stationary solutions ^1,2 (a;) are 
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FIG. 11: Phase diagram of refraction angles as a function of V on x-axis and (p on y-axis, and 
for an MVP barrier of width D = lOOnm at B = 0.1 and 3T. The gap corresponds to the region 
where TIR occurs. 



^2 



(30) 



These are similar in form to those for a pure magnetic barrier or an electrostatic step 
But [kx^ ky) and [qx^ky + j-]^ namely the x and y components of the wavevec- 



potential 



10 



tor inside and outside the barrier regime are different. Here, s, s' are sgn{E — V) in the 



respective regions. Upon setting incident energy Ep as hvpkp, substitution of Eqns.(29) 



and (30) in Eq.(28) leads to 



k^ -\- k 



ql + {ky + 



p, with kp 
V 



Ep^ 

hvi 



\x\ > d 



{kp- 



k p ^ I X I ^ d 



(31) 



Ib hvp 

The incidence angle and the refraction angle 6 are given by tan-i(f2^) and iBjr\^-^^) 



respectively. Eq.(31) can then be rewritten to obtain the Snell's law analogue for electron 
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waves of such massless Dirac fermions incident on the EMVP barrier as 

sin 1^1 = ^i.(sin|0| +sgn(0)-^) (32) 



Sp sin \6\ 



Sp = ^ = ^l^^l; s =sgn{Ep-V) 



Comparison between Eq.(23) and Eq.(32) shows that the potential barrier effectively 



scales the refraction angle by the scale factor Sp defined above. Sp is a non-monotonic and 
discontinuous function of V [FigjToJ^a)] and we shall study its impact on the refraction of the 
incident electron wave. There are basically two regimes to be analysed depending upon the 
strength of electrostatic potential applied (i) V < Ep, (ii) V > Ep. For positive incidence 
angles = (0,7r/2), in the absence of any electrostatic potential the refraction angle \6\ is 
larger than the incidence angle |0|, thus the electrons are seen as passing from denser to rarer 
medium. As V increases from to Ep, the function 5*^ and hence 6 increases continuously 
(this is shown in figjlOta) and (b) respectively) and eventually the electron wave suffers 



TIR as soon as the r.h.s of Eq. 32 becomes greater than 1. Thus, by increasing V it is 
possible to totally reflect an electron wave for any given and B. Since the reflectivity of a 
magnetic barrier increases with higher B, this implies that the addition of V can effectively 
convert a weaker magnetic barrier into a stronger one. This effective control of the strength 
of the magnetic barrier by an electrostatic potential ( gate voltage) is a consequence of the 
ultrarelativistic nature of the charge carriers in graphene. In a subsequent section we shall 
provide a more formal reasoning for this feature. 

When V surpasses Ep, the sign of sin \ 6\ will be opposite to the sign of sin |0|. At smaller 
V close to Ep there is still TIR, whereas for very high V much above Ep the refraction 
becomes negative and the wave again retraces its path back in the barrier regime. This 
situation is depicted in Fig. [Io|(b), where the wavevector for the refracted ray is shown 
changing with increasing V. A similar analysis can be done for negative incidence angles. 
We summarize this discussion by plotting sin \ 6\ as a function of the incidence angle for 



different V and B values in Fig. 11 As can be seen, each of FigjTTj^a) and (b) are separated 
into an upper and a lower part by the line (ps = ^^^~^{~'k^)- In the lower part where |0| is 
larger than |0s|, with increasing V the absolute value of the refracted angle, \d\ exhibits four 
phase regions: rarer — )• denser, denser — )• rarer, TIR gap, negative refraction. In the upper 
part, the four phase regions exhibit a different order with increasing V: negative refraction, 
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FIG. 12: Here and ky represent the transverse and longitudinal components of the wave vector. 
The circles represent solutions of the dispersion relations in the domains 1 and 2 and correspond 
to propagating waves. [82] 

TIR gap, denser — )■ rarer, rarer — )■ denser. All these different regions meet at the limiting 
point {EF,(f)s), where the behaviour is singular and is discussed later on in this article in 
more detail. 

The refraction laws for the electron waves incident on an MVP/EMVP barrier can be 
alternatively understood with the following geometric consideration |82] . We consider two 
domains with scalar and vector potential {Vi,Ai) and (^2,^2) respectively. The electron 
wave being incident from domain 1 at an angle and refracted in domain 2 at angle 6. Then 
the parallel and perpendicular wavevector components in domain 1 and 2 are related by the 
following dimensionless dispersion relations: 



'^1,2) 



(33) 



Each of the eq. 33 represents a circle in the (k^^ky) plane, so that the wavevector for 
the propagating solutions will lie on the circle itself. Now, the refraction laws can be easily 
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derived from fig. 12 as follows. Due to translational invariance along the y-direction, the 
perpendicular wavevector component takes the same value in both domains ky-^ = ky^ = 
kyg. For a given ky^, propagating solution kx in domain 1 and 2 will take the value which 
lie at the intersection of line ky^ = constant with the corresponding circles. There are 
two such intersections for each of the circle. (This is shown in figure at Oi^2 and 0[2)- 
The physically meaningful intersection will be determined by the sign of (e — 'Ui,2)- For 
definiteness we assume the incident current density to be positive i.e. k^^ > 0. Then 
using the parametrization k^^ = (e — Mi)cos0, cos0 will be positive or negative when 
{e — Ui) > or < respectively. This means the incident wave will be directed along the 
vector joining the centre of circle with Oi when {e — ui) > or opposite to the vector 
joining the centre of circle with O'l when (e — ui) < 0. Similarly, using the parametrization 
kx2 = {^ — ^2) cos 6, the refracted current (which should possess the same sign as the incident 
current density) will be directed along the radius vector from the centre of the circle to the 
intersection point O2 when (e — U2) > 0, or opposite to the radius vector from the centre 



of the circle to O2 when (e — U2) < 0. Simple geometry from the fig. 12 gives the following 



relation that defines the connection between the incident (p and refracted 6 angles 
kyg = —Ai — (e — Ml) sin = —A2 + (e — U2) sin 9 



This relation gives the same condition as in eq. 32 for EMVP barrier. 



2. EMVP barrier: Transmission probability 

Transmission through a magnetic barrier gets strongly affected due to the presence of V. 



Again from the continuity at x = ±d in Eqs.(29) and (30), the transmission coefficient can 
be obtained as 

2ss'e~*''^'^ cos0cos^ 
~ gg/jg-ig^D (.Qs(^0 + 9) + et^^ cos(0 - 9)] - 2i sin q^D ' ^ ^ 

where D = 2d and = k'p cos 9. Here s and s' corresponds to sgn(E) and sgn(E-V) 
respectively. The expression for t is same as the one for electrostatic potential barrier [TT] or 
MVP barrier. However the difference in the result is due to the change in the expression for 
q^, namely the x-component of the wave vector inside the barrier regime. In figjT3]is plotted 
transmittance T = t*t as a function of and V in the presence of an EMVP barrier. When 
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FIG. 13: Transmission in the presence of EMVP barrier in monolayer graphene, B = IT, Ep = 
50mey, D = lOOnm. 

the potential is lesser than but close to Fermi energy, transmission takes place over a very 
small window located asymmetrically along the axis. As the point V = Ep is crossed, the 
function Sp changes sign thereby changing the sign of refraction angle. Thus, the potential 
barrier induces more asymmetry in the transmission as against pure MVP barrier, this is 



clearly depicted in Fig. 13 



The point V = Ep represents a singularity in the spectrum and demands a separate 
discussion. In the absence of magnetic barriers, such a point represents the zero modes for 
Dirac operators and leads to the emergence of new Dirac points. This has been discussed in 
a number of recent works considering Andreev Reflection in a graphene based NIS junction 
[S3] and for electronic states of graphene in a periodic potential [S51 - I57] . The presence of a 
magnetic barrier breaks the time reversal symmetry explicitly and these zero modes become 
the zero modes of the modified Dirac operator and the corresponding solutions are different. 
The equations satisfied by ■ipi^2 are 

[d, - {ky + 1/Ib)] = 

[d,, + iky + l/lB)]ij2 = 

for which the solutions are obtained as 

^2 = c2e-('=^+^/'^)^ 
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FIG. 14: Averaged transmission Vs voltage through a single EMVP barrier in monolayer graphene 
for different magnetic fields. 

By using the continuity condition at x = ±d, the transmission expression is obtained to be 



The significance of such zero modes is that on the two sides of the singular point, the relative 
sign between the and 9 becomes opposite. 

3. Conductance: Effect of various EMVP Barriers on Transport 

To see how the above angle-dependent transmission properties effect electron transport 
through such barriers, we plott the average transmission through the barrier as a function 
of the potential V at various strengths of the magnetic barrier. The average transmission 
at a given barrier strength V and B is defined as [73] 



This formula, when generalized to a range of energy levels, leads to the Landauer con- 
ductance G oc 2TTne'^{T{B,V)) /h. This has been plotted for three different strengths of 
magnetic barrier B over a range of potential barrier strength V for a single EMVP barrier 



t = 



2cos0e ^''^ 



(3i't>Q{ky + l/lB)D _|_ Q-i4>Q—{ky+l/lB)D 




(35) 
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in Fig. 14 The transmission shows a minimum as expected when the point qV = Ep is 
approached. The transmission grows on both side of this singular point and finally at higher 
\V\ oscillates around an average value. 



4-. Generalisation to MVP/EMVP harriers of any arbitrary shape 

A general algorithm for the calculation of electron transmission in graphene through 
inhomogeneous electric and magnetic fields of any abitrary profile has been carried out in 
[88] using transfer matrix method. The fields are invariant in one direction; and the method 
involves the division of the one-dimensional domain into slices and taking an appropriate 
approximation of potential form in each slice. The equation for each slice is solved and the 
continuity conditions are used at the interfaces of two such slices. The exact solution of the 
equation in each slice depends on the potential form chosen. It is shown that in the presence 
of piecewise constant scalar potential and piecewise linear vector potential the resulting 
equation admits parabolic cylindrical functions as the solution basis in moderate magnetic 
fields while the basis functions tend to complex exponentials in the presence of extremely 
small magnetic fields. Also it is shown that the localised charge distribution (electron and 
hole puddles) which arises due to the presence of disorder in graphene [TS] can be modelled 
as the charge as arising from a scalar potential distribution proportional to it. With this the 
transmission calculation corresponding to a one- dimensional potential extracted from the 
experimental data in ref. [78j is presented. The transport analysis with magnetic barriers 
in which the edges are smoothed out have also been carried out |89]- where a hyperbolic 
profile has been chosen. Here the corresponding Dirac equation can be analysed within the 
formalism of supersymmetric quantum mechanics, and leads to an exactly solvable model. 
Also, exact solutions have been given for a Dirac electron in the presence of an exponentially 
decaying magnetic field [90] • 

C. Quantum Goos-Hanchen Shift in single MVP and EM VP barriers 

Analogues of optical phenomena such as refraction [TT], coUimation [91], Fabry Perot 
interference ^2], Bragg refiection [19] etc. in the ballistic transport regime for Dirac fermions 
in mono and bilayer graphene based structures in presence of scalar and vector potentials 



35 



have been proposed. We discuss them in various sections of this review. In the following 
section we discuss another important optical phenomenon, the Goos-Hanchen (GH) effect 
[nSl [M] • The GH effect is a phenomenon of classical optics that describes the lateral shift 
between the centre of a reflected beam and that of incident beam when a total reflection 
occurs at the interface between two media. The shift occurs as the totally reflected ray 
undergoes a phase shift with respect to the incident beam, and this is is detectable since the 
extent of a real beam is always finite. The shift reverses sign if the second medium behaves 
like a metamaterial with negative refraction [^1^. Such a lateral shift for totally as well as 
partially reflected electron waves can also occur for non-relativistic electrons passing through 
a semiconductor barrier [97], magneto-electric semiconductor nanostructure [98j. Recently, 
it has been shown that ballistic electrons passing through a p-n interface in graphene [99] 
also suffer a GH shift, which changes sign at certain angles of incidence. The analysis was 
extended for the case of MVP and EMVP barriers subsequently [21]. To understand GH 
shift in graphene, we consider the following wavepacket (beam) of electrons impinging on a 
p-n interface or a MVP /EMVP barrier at energy Ep: 



dkyfiky - k)e'^yy+'''^(^y)^ 



(36) 



The envelope function ensures the wavepacket is of finite size along the y-direction and is 
sharply peaked at ky = k. Thus, k G {Ojkp) and the angle of incidence (piky) E (0, |). 
This fact is represented by writing the x-component of wavevector, k^ as well as (p both as 



function of ky in Eq.(36). We take a finite beam with a gaussian envelope such that 



fiK -k)= exp[- 



2A^ 



(37) 



When Afc ^ kp, we can approximate the fcj^-dependent terms by a Taylor expansion around 
k and retaining only the first order term to get 



(j){ky) « (j)ik) + gj^lkiK - k); k^{ky) ^ k^{k) + -Q^VkiK - k) 



(38) 



Substituting in Eq.(36) and integrating, we obtain 



Hy-yT? 



e 2 



(39) 
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where 

= -KCk)x, r = -KCk)x - <P'Ck) (40) 

Here the primed (/) quantities denote derivative taken wrt ky. Thus, upper and lower 
components of the spinorial wave function are locahzed at separate points along the y-axis. 

The reflected wavepacket can also be written in an analogous way by making the trans- 
formation kx to —k^ and to vr — as well as multiplying the reflection amplitude 

r{ky) = \r{ky)\e^'^''^''y\ The reflected wave is then 



dkyfiky - k)e'^yy-'^-^^y^''r{ky 



(41) 



Here again, s = sgn{Ep). The spatial profile of the reflected wave can be again obtained 
by first expanding all ky dependent quantities around k and retaining only the first order 



terms and then integrating in Eq.(41). This leads to 



= ^J2'K^le''^^y-^^^~^^\{k)\ 



e 2 



—se 2 



(42) 



Here, y\. and y'[_ are given by 

f+ = -<P'r{k) + k'^{k)x,f_ 



-<i)',{k) + k',ik)x + (P'{k) 



(43) 



The above expression shows that the upper as well as lower components get shifted because 
of the phase factor. The GH shifts of the upper and lower components are respectively given 
by 



a+ = yl-y^ = -<P\ik) + 2k'Mx 



y- - y- 



h'{k) - (t)',{k) + 2k',{k)x 



Thus, the average shift for an MVP or EMVP barrier is 

(T = ^(a+ + a_) = ct>'{k) - 0;(fc) + 2k',{k)x 



(44) 



(45) 



The situation is depicted schematically in Fig. 15 The last term in the above expression is 
a coordinate dependent quantity and will get an equal and opposite contribution from the 
—(f)'^{k) term. The resultant a will thus be independent of the choice of the coordinate of 
the interface from which TIR will take place. 
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FIG. 15: GH shift for an EMVP barrier. The solid hne corresponds to the upper component of the 
pseudospinor and the dotted hne to the lower component. The GH shift can be either (a) positive 
or (b) negative. 

1. GH shift at p-n interface 

To calculate the GH shift occuring at the p-n interface, we first require to calculate the 
gradient of the phase of the refletion coefficient, the negative of which gives the GH shift. 
For this we match the wavefunctions at the interface (x = —d) of the barrier such that 
the propagating wavefunction on the left of the interface should match with the evanescent 
wavef unction at the right of the interface. 



^r = <^ (46) 



e 



_|_ ik^x 2; < — C? 



a'e 



-K{x+d) if x> -d 



gl^i{k^X + <P) _ ^^-i{k^X + <P)l if X<-d 

^2 ={ (47) 

-■i7s'a'e"''(^+'^) if x > -d 



Here s = sgn{Ep), s' = sgn{Ep — V), and 



hvpiK + ky) r, [E — V 

7 = C — 

' T^. — M ' y 



2 



E-v ' y \ hvF 

TIR at p-n interface will take place when = E'l' = {Ep — V^)^ > 0, for which the incidence 
angle lies in the range (p > (pc = sin~^ — 1^ The continuity of the wave function at 
X = —d gives the reflection coefficient as 

r = e-^-^[ ] = exp(-zfc.D) exp(2i^) 

with tan 5 = tan + ss'7 sec (48) 
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FIG. 16: Dependence on the angle of incidence (p of the GH shift a, for Uq/Ep = 1.5 (dashed 
curve, p-n interface) and for Uo/Ep = 0.4 (dotted curve, n-n interface). 



Clearly, the phase of the reflection coefficient is (j), 
written as: 

a = (jy'ik) — 2(5'(fc), where 5 = tan~^(tan0 + ss'7sec0 
Using this expression, one can calculate the GH shift at p-n interface. 



—kxD + 25. The GH shift can now be 

(49) 



, . (50) 

n sm (j) cos (p 

This is the same expression as the one obtained in [99] . Since /t > for decaying evanescent 
solutions at the right of the p-n interface, for positive incidence angles the negative GH shift 
at the p-n interface {Ep < V) will occur when sin^ — (V/Ep — 1) < 0. Thus, for angles of 
incidence 0c < < 0* = sm^^v^sm^^, GH shift occurs in the backward direction. This is 



explicitly shown in fig, 16 



2. GH shift for MVP/EMVP harriers 

For pure MVP barrier, TIR will take place when \sin \9\\ = sin\(f)\ + sgn{(f))-^^ > 1. 
This will happen when sin |^^| > 1 for < < | and when sin |^^| < —1 for — | < < 0. In 
the latter case, this requires the wave vector to be negatively refracted at sufficiently high 
magnetic field before TIR occurs. This is already explained in section IV A in detail. The 



GH shift can then be obtained with the expression as in eq,49 with only 



7 



K,+ {ky+ 
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Ep 

hvf 
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FIG. 17: (color online) GH shift for an EMVP barrier at (a) O.IT and (b) 3T. The y-axis corre- 
sponds to the incident angle (p. The green regions have propagating solutions with no shift. 

An EMVP barrier is essentially a p-n junction placed in a MVP barrier. For an EMVP 
barrier, TIR occurs for electrons incident from both sides of the surface normal but at 
different critical angles. For a given 0, it is possible to change V adiabatically and get the 
electron wave reflected over a range of V satisfying \sin \ 6\\ > 1. We can keep (p fixed and 
increase V. At V = V^^, TIR occurs when sin |^^| = 1 giving 

^ = 1 - [sin + ■'^9n{(l>)-^^] 

Upon further increasing V, the electron wave remains totally reflected till V reaches the 
second critical value Vc^ such that sin \ 6\ = —1 giving 

— 1 + [sin 101 + sgn{(f))- '' 



F^B 



TIR occurs in the range V G [Fci, ^^2]- Fig. 17 is plotted the GH shift over the entire 
range of V and 0. At all other regions in the ip — V plane the GH shift is set to (green). 
The explicit expression for the GH shift for an EMVP barrier in dimensionless form is 

1 -2ff cos2(5(l + tan0tan5 + ss'2^) 
afcp = '-^ (51) 

COS0 

The ^23Y term in the numerator which is coth a or sec /3 diverges at the critical angle for TIR 
and leads to a divergent GH shift. This is because just at the critical angle the wavevector 



lies in the interface of the two regions. As a result in Fig. 17, the border TIR region of finite 
GH shift shows the highest GH shift. 
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The upper and lower parts of the curve are again separated by the hne 0^ = sin~^ (^fcAe) 
. The left and right boundaries correspond to < Ep and V > Ep- In the lower part, the 
GH shift is mostly negative with the left and right boundaries having positive and negative 
refraction respectively. 



Comparing Eq.(51) with Eq. (50), we see that 7^ — cr(— 0). This is due to the non- 
specular nature of electron refraction at an EMVP barrier. This is an important difference 
in the quantum GH effect that occurs upon TIR by an EMVP barrier as compared to TIR 
by a purely electrostatic barrier for which \a{ 



|cr(— 0)1 as can be seen from Eq,50 



3. Lateral shift in the transmitted wave 



As explained in the previous section, Goos-Hanchen shift is analogous to the phenomenon 
of the lateral shift of the light beam total internally reflected from dielectric surface. Chen 
et al. [971 [98] investigated the lateral shifts for Dirac fermions in transmission through 
monolayer graphene barrier, based on tunable transmission gap. This shift has same physical 
origin that is due to the beam reshaping since each plane wave component undergoes different 
phase shift, however it has nothing to do with the evanescent waves which plays an all 
important role in the lateral shift of total internally reflected wavefunction. For this reason 
this shift can be termed as Goos-Hanchen like (GHL) shift and can be considered as an 
electronic analogue of the lateral shifts of light beam transmitted through a metamaterial 
slab. 

According to stationary phase approximation |100] and on the same lines as explained 
for GH shift for the reflected beam in th previous section, the GHL shifts of the transmitted 



beam through a barrier in graphene can be obtained as [HZl [ 



crt = - 



Sky 



where (pt refers 



to the phase of the transmission amplitude and the subscript corresponds to the wavevector 



at central incidence angle. Using (11), and assuming the form t = Itle*"^', the transmission 
phase 0t is obtained as 



tan 



ss sec 9 sec — tan 9 tan ) tan q^d 



so that 



crt 



5^ 
6L, 



(itan0n ■ T 



, where /cg = ky^ — ss'kpk'p. 



^2 ^2 



sin Iq^d k\ 
2q^d ql 



(52) 
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FIG. 18: Transmission gap (a) and GHL Shift (b) as a function of incident energy, E at 
U = 62meV, d = lOOnm. 



30°, 



Since the product ss appears in the GHL shift expression so the analysis of the expression 
can be classified into ss' = —1 (Klein tunneling regime) and ss' = +1 (classical motion 
regime) . 

In the regime when ss' = — 1, /cq = kpk'p + k^^ > (for (f) > (f)c and for all values of 



incident energy, E). Then in accordance with (52), the maximum absolute value GHL shift 
corresponds to Qxd = nir (These wavevector also corresponds to transmission resonances). 
Hence 



n 7 

2~"t^^' 

1x0 



Thus in Klein tunneling regime, the GHL shifts are negative shifts when </> < (p^. 

In the other regime when ss' = 1, /cg = fc^g — ss'kpk'p = —kl^ < (for (j) < (j)c and for all 



values of incident energy, E). Then in accordance with (52), the maximum absolute value 
of GHL shift corresponds to qxd = rni (This wavevector also corresponds to transmission 
resonances). Hence 



q^d=mT 



k'^ 
1x0 



(itan< 
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FIG. 19: (Color online) (a) Contour plot of the transmission through a single barrier Vs its height 
V and wave vector ky for D=100 nm. (a) B = (b) B = O.IT. 

Thus in classical regime, the GHL shifts are always positive. 

Thus we see that the GHL shifts in transmission can be negative or positive, and can be 
enhanced by the transmission resonances when the incidence angle is less than the critical 



angle for total reflection. This is clearly depicted in Fig, 18 



D. Fabry-Perot Interference in Graphene Heterojunctions: Effect of magnetic 
field 

Another important optical analogue that can be studied in the ballistic transport regime 
of graphene is the formation of Fabry-Perot resonances. Study of such resonances is the linear 
and non-linear transport regime for non relativistic electrons for devices such as resonant 
tunnelling diodes is well known for a long time particularly after the pioneering work by 
[H 1101] . A excellent review on the experimental progress on such resonances in graphene 
based heterostructures in the ballistic transport regime can be found in [15] . 
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Below we revisit the analysis of potential barrier as a Fabry- Perot interferometer [HI [92] 
and then we analyse that how the FP resonances get effected in the presence of magnetic 
field in graphene. 

The two edges of the square potential barrier acts as the two interfaces of FP interfer- 



ometer p2\ , as shown in Fig. 20 As an electron wave is incident on the potential barrier at 
an angle 0, it splits into the transmitted wave and reflected wave at the left edge {x = —d). 
The transmitted wave propagate further to suffer another refiection(transmission) at the 
right edge (x = +d). In this way it gets multiply reflected between the two edges of the 
barrier. The difference in the optical paths between the transmitted wave ti and t2 which 
suffer consecutive reflections at the right interface is given as 

AL = 2(1 - V/E){BC + CD) - BP 

V\ 2D 

1 - 2D tan 9 sin 6 

E J cos 9 ^ 

^ V\ 2D 2D y\ ■ 2a 

1 1 sm'' 9 

E J cos9 cos9 V EJ 

V\ ( V\ 

1 — — j 2D cos 9 using sin = ( -'- ~ ^ j ^ 

This means that the net phase difference between ti and ti'vs b = k^^L = 2qj.D. 

Again we introduce r and t to be the reflection and transmission coefficient for the 
potential step outside the barrier, and r' and t' the corresponding coefficients for inside the 
barrier. In analogy with the optical waves, net transmission through the barrier is obtained 

as 

ttot = tt' + ttV'e^^ + . . . + ttV''^"-')e*("-i)'^ = tt'/\l - r'\'^] 



Then, the total transmission probability T = t^^t^tot is obtained as 



T = 1/ 



1 + 



4|r 



/|2 



'1 



r'\2\2 



sm^6/2 



(53) 



Proceeding with the calculation given in section 2, the reflection and transmission coef- 
ficients r, r', t and t' can be obtained by matching the wavefunction at the interface of the 
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potential step. 

V'x<0 = 

i>x>o = 

Matching conditions at x = gives 

t = t' 



2 COS I 



1 — s'/s cos(<; 
1 + s' /s cos(<; 



(54) 



Substituting [54] in [53] we obtain the total transmission probability of a single barrier. 

^ cos^ 9 cos^ (b 

T = - 

{cos q^D cos 9 cos 0)^ + sin^ QxD^l — ss' sin 9 sin0)2 

Clearly the transmission resonances which are also the Fabry perot resonances occur at 
QxD = nn. The energies at which these resonance occur are obtained as: 



E = V±hvF[n\yD^ + eW^ 



(55) 



Also T is equal to 1 at normal incidence = 0. 



To show the effect of magnetic field which is already mentioned in section II B explicitly, 
we consider a potential barrier of height V with a commensurate _L magnetic field B [92] 



V{x) = V, B = B&{£ - x^)z \x\ < d 



(56) 



With the similar analysis as in section ]III B[ the carriers in the presence of such potential 
form, obey the Schrodinger like equation of the form: 



92 



{kyls + X 



-\2 



(e — f )^ ipi^2- \x\ < d 



1,2 



Upon solving, the solutions can be obtained as: 



{V2{x + kyle)) + 6L'(e-i;)2/2_i(-y2(x + kyle)) 
^2(3;) = a^^— ^£'(e-{))2/2(V2(^ + kyls)) + b-^—-^D^^_i)2/2{-V2{x + kyle)) 
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B=0 B^O 

FIG. 20: Electron transmitting through a potential barrier in (a) absence (b) presence of magnetic 
field. 



The solutions in the outside regions namely |x| > remain same as in section IIIB By 



using the continuity condition to match the wavefunctions at the boundaries of the barrier we 



numerically evaluate the transmission as plotted in fig, 19 Since in the presence of magnetic 
field the y-component of the momentum varies inside the barrier as ky ky — eBx/h, so 
the for the incidence angles at the two interfaces 1 and 2 to be of equal sign, ky{xi) and 
ky{x2) should have opposite signs ie. 

ky — eBxi/h = ky + eBd/h < 
ky — eBx2/h = ky — eB{d)/h < 

This gives the following condition on the y-component of momentum: 

- eBd/h < ky < eBd/h 

A fringe contrast from color minimum of transmission to the color maximum of the trans- 



mission on the two sides of the lines at ky = ±eBd can be clearly seen in the Fig, 19 which 
occurs as a consequence of the half a period shift in the transmission resonances. Occurance 
of such phase shift for charge carriers in graphene was originally studied for a harmonic 
potential in a uniform magnetic field by Shytov et al. |H] and used as an experimental 



signature of Klein tunnelling [3l] . This has already been discussed in section II B Trans- 
mission and conductance of massless Dirac fermions in the presence of multiple, disordered 
short range scatterers can also be understood in terms of resonant transport through double 
barrier structure. This has been discussed in detail in Ref. 
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V. PERIODIC LATTICE OF MVP BARRIERS 



Electron transport in graphene in presence of superlattices formed with electrostatic 
potential or magnetic barriers lead to large a body work in recent times UTj I1U2H104] . 
They lead to interesting physics and potent with the possibility of creating new devices in the 
ballistic transport regime since they significantly alter the band structure. On other hand 
they resemble closely the phenomena of light propagation through medium with periodically 
modulated refractive index/dielectric constant which has been studied in detail in optics 
literature and also led to many important applications [31 1105] . In the following section we 
discuss this type of transport and their consequences. We shall also discuss the emergence 
of additional dirac points in these modified band structure and the resulting coUimation of 
transport electrons using such superlattice structure. 

A. Bandstructure modification 

In the previous section we have shown how MVP/EMVP barrier affects the ballistic 
transport in a profound way. Thus it will be interesting to find out how a periodic arrange- 
ment of these barriers will effect the transport, what is the resulting band structure etc.. As 
mentioned earlier, it has a direct optical analogue, namely with the electromagnetic propa- 
gation in periodic stratified media following the pioneering work by Yariv and collaborators 
[3] . It also provides us relativistic version of ordinary Kronig-Penny model and its magnetic 
counterpart. 

Practically there will only be a finite number of barriers present in the system and the 
lattice translational symmetry will break down at the boundaries. However to simplify 
the analysis, we assume that the unit structure- a double EMVP barrier structure can 
be repeated infinitely along the x-axis. The magnetic field creating such a barrier B = 
Bz{x)z = BIb[S{x + d) + 6{x — d) — 2B5{x)]z together with a superimposed splitgate voltage 
geometry, hence forth will be called DEMVP barrier. Since changing or reversing V can 
locally convert a charge-neutral region into a p-n or n-p junction, such a periodic barrier 
can also be thought of as a semiconductor heterostructure. It should be emphasized that 
this can lead to new device structures due to combined effect of the highly inhomogeneous 
and periodic magnetic fields and controllable voltages. The vector and scalar potentials are 
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characterized by 



V{x) = -V-Ay{x) 
V{x) = V-Ay{x) = 



?B, —d < X < : region I 
?B, < a; < (i : region II 



(57) 



The schematic diagram in Fig. [9[b) shows how such DEMVP barriers can be constructed 
by placing metal electrodes and insulating layers between PMA materials with alternating 
magnetizations and the graphene sheet. 

We consider each unit cell of size D = 2d for the MVP as well as for the electrostatic 
potential barriers. Thus, the n-th cell is given by the region between (n — 1)D and nD. In 
the a-th part of a unit cell, the wavefunction is 



(58) 



Here, 



a 



Sn = Six, si,2 = sgn{E ± V) q'l 2x = 9i,2 



(59) 



The wavevectors q{i,2) are given as 



ql + [ky + 
ql + [ky 



1 , 



Ib hvp 

The exponential factor e~"^ reveals the existence of lattice translational symmetry, which 
is not present for the isolated EMVP and DEMVP barriers. 

The continuity of the wavefunction at the first interface aX. x = {n — 1)D gives 



1 



1 



Cn-1 
dn-1 



i[qiD-e 





an 




bn 



(60) 



Similarly, the continuity at the second interface ai x = [n — 1)D + d gives 







an 














bn 








dn 



(61) 
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(a) B = O.ir, gV/£f =-0.7 
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a = Q.ir, ijl'/Bi- = 0.7 
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FIG. 21: Band structure at some characteristic values of V with i3-field strength strengths (a) 3 
Tesla, and (b) 0.1 Tesla. 



Introducing a shorthand notations like, 



Ma 



M. 



e*'^- 

1 

S, Si2 



Mb,, 



Me,. 



1 1 



Now Eqns.(60) and (61) can be rewritten as 



Ms.„.,Mg, 



Ms,Me,MB, 



Cn-l 
dn-l 

an 
hn 



Ms,^Me,MB: 



Ms,Me,MB, 



an 

K 

dn. 



The above two matrix equations can be combined as 



(62) 



' M^^'M-'Ms,Mg,MB,M^^'M;^'M,2Me,MB, 

dn-l 

According to Bloch theorem, 



Cn 




Ku 






Cn 


dn 




K21 


K22 




dn 



Cn-1 




Cn 








dn-l 




dn 



(63) 



(64) 
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The matrix Kmat = 
eigenvalue equation 



K21 K22 



is unimodular. From Eqs.(63) and (64) we obtain the 





Kl2 








Cn-l 


K21 


K22 




dn 




dn-1 



(65) 



where K is the Bloch momentum. The complex conjugate eigenvalues A are given by 
det\Kmat — A/| = O.This implies Ai + A2 = exp{—iKD) + exp{iKD), which finally gives 



Ki^,B) = ^cos''[^TriK„,at)] 



(66) 



Writing in terms of the wavevectors qi,q2 and the angles 61,62, the above eigenvalue 
condition reads 



cos KD = cos qidcos q2d + sin qid sin q2d x 



tan 61 tan 62 



S1S2 



(67) 



cos 6*1 cos 62 

This equation provides the band structure for a periodic DEMVP barrier superlattice in 
general for any V E. By substituting V = in the above expression, the band structure 
for a periodic DMVP barrier superlattice can be obtained. This is an extension of the Kronig- 
Penney (KP) model to two-dimensional massless Dirac fermions. Thus, it is interesting to 
compare the DEMVP band structure with other variants of the KP model. The original KP 
model describes Bloch waves in a one- dimensional periodic potential Several authors 
have also studied the relativistic KP model |106H109] , where the motion considered is strictly 
one-dimensional. The non-relativistic KP model in periodic structures created by MVP 
barriers has also been studied jlSI UlllOj . For graphene charge carriers, such problems have also 
been studied for different types of periodic magnetic [191 El], 11031 1104j as well as electrostatic 
[102] barriers. Magnetotransport studies in presence of periodic barriers |lllj where the 
periodic modulation is sinusoidal was carried out for monolayer graphene and it was found 
that such periodic modulations Landau levels into bands which oscillate with the strength 
of the magnetic field. 

In Fig. [21] is plotted the band structure of corresponding to the periodic DEMVP barriers 
as a function of magnetic field B and voltage Conducting regions are seen over a wide 
range of E with forbidden regions in between. At O.IT, the region near ^ = is mostly 
conducting for different values of energy where a forbidden region starts opening up both to 
the left and right of the V axis. The gaps are much wider than their counterparts for pure 
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FIG. 22: The group velocity Vg (in units of vp) at the neutraUty point as a function of d/lB- 

MVP barriers. Close to qV — ?■ Ep, an extended forbidden zone appears. To understand 
more explicitly, the band structure is also plotted at ^ =0.7 close to the singular point. 

Here, the conducting regions are intervened by large patches of forbidden zones. In 
comparison, at = —0.7, the behaviour is completely changed and the forbidden zone over 
the same range of kyls shrinks considerably (left column). For higher V where qV > Ep, 
the system is conducting at almost all E. At 3T and near zero V, a forbidden region opens 
up at various values of the energy and is much larger than at O.IT. The gapped regions 
to the left and right of zero V at 3T are located in a pronounced asymmetric manner as 
compared to when B = O.IT. This asymmetry in the band structure as a function of V 
as well as the opening up of large forbidden zones for certain values of V differentiates the 
transport through EMVP barriers from MVP barriers and thus, it provides more flexibility 
to tune such transport. 

As shown by various studies, carriers in graphene superlattices exhibit several interesting 
features [112] |113j [91] that result from the particular electronic SL band structure. Below 
we briefly analyse the bandstructure of periodic DMVP barrier superlattice. 



1. Dispersion in the vicinity of neutrality point 



We begin with analysing the dispersion relation in the vicinity of the neutrality point i.e.. 
close to zero energy. By rewriting the Eq,67 in the form of f{E, k^, ky) = 0, then expanding 
the /-function to the lowest order in E and ky and then in k^ in the vicinity of the neutrality 
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FIG. 23: The circles emphasize the finite energy Dirac points at the values corresponding to Eq. 
691 for S = O.Sr, and d = lOOnm. 



point {E = 0, ky = kx = 0), we obtain the dispersion relation: 



d/l 



B 



E = ±hvFVox kl + k'^, with vo = . 

^ sinh(a//_B) 



(68) 



Clearly, the above dispersion relation is isotropic {dE/dk^ = dE/dky) with zero gap between 
the valence band and conduction band. Moreover the group velocity Vg = VqVf depends on 
the width of the unit DMVP barrier as well as on the 5— field strength. The group velocity 



is plotted in Fig. 22, which clearly shows that Vg is always smaller than the Fermi velocity 



and it monotonously decreases for increasing d/lB-, which is implied from the single MVP 
barrier analysis as well. 



2. finite energy Dirac points 



Apart from the neutrality point described in Eq. 68 , many other degeneracy points can 



be seen in Fig, 23 where the dispersion have a double cone like structure. The locations of 
these points occuring along the ky = direction can be obtained analytically as 

En = ±hvF^{lllBY + {mi/dY (69) 
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FIG. 24: Plots of the dispersion relation at n = 1 (top row) and n = 3 (bottom row) for B = O.IT 
(left column) and B = IT (right column), d = 50nm. 

Again by expnading the /-function in the vicinity of these points En{kx = 0,ky = 0), the 
dispersion relation is obtained as 



E = E^±hvF^vl,kl + vlykl n = l,2,... (70) 
with the group velocity components given by 



fn7r)2 



2 



{d/hy + (nvr) 
{d/h? + (nvr)2 

Clearly the dispersion relation is linear in k^ and fcy, which shows that in the presence 
of DMVP barrier superlattice, there are infinite number of degeneracy points close to which 
the dispersion has a double cone shape which is asymmetric in k^ and ky direction. These 
points are termed as "Extra Dirac points" in literature |1 14] PT5] . 

In contrast to the dispersion relation at the neutrality point which is isotropic in nature, 
the dispersion at these finite energy extra Dirac points is anisotropic [dE/dk^ 7^ dE/dky). 
As determined by EqjTO} we plot in Fig, 24, the dispersion relation at n = 1,3 for different 
5— field strengths, B = 0.1, IT. From this we infer a collimation along ky direction ie. 
Vy oc dE/dky vp and f ^ as the contours become almost independent of k^. This is 
similar to that found for a superlattice of electrostatic potential barriers 
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VI. THE TRANSFORMATION OF MAGNETIC TO ELECTRIC FIELD FOR 
DIRAC FERMIONS IN GRAPHENE: A CASE FOR RELATIVISTIC INVARI- 
ANCE 



For the electron transport in presence of a transverse uniform magnetic field it was shown 
by Lukose et al. |116] that the Landau levels gets modified when the system is subjected to 
an electric field, finally leading to the Landau level collapse through the renormalization of 
the cyclotron frequency. A similar rescaling of the magnetic field strength took place when 
a suitable gate voltage is applied for the magnetic barrier problem as well [21] as discussed 



in section |IVB[ To summarize one can change the strength of a magnetic field by suitably 
applying a get voltage. A more detailed analysis done later on |117] how the relativistic 
invariance of the effective equation that governs the charge carriers in graphene leads to 
this interconvertibility of electric and magnetic field. In this section we shall describe this 
approach. 

According to Special theory of Relativity, the laws of physics must be formed in a form 
so that it do not allow one to distinguish between the frames of reference which are moving 
with a constant relative velocity. The transformation between such frames are described by 
Lorentz transfromations. Like Maxwell's equations, Dirac equation is also Lorentz covariant. 

So, for the charge carriers in graphene obeying "Dirac-like" equation, the Lorentz covari- 
ant structure of Hamlitonian can be used to find an analogy with special relativity. The 
analogy extends to problems involving coupling of electrons to external electric or magnetic 
fields. It is known from from special Relativity that in the presence of electromagnetic fields 
there are two invariants; namely E'^ — and E.B which remain unchanged in transition 
from one inertial frame to another |118] . This means, by means of Lorentz transformation 
we can always give E and B any arbitrary values, subject to only condition that E'^ — 
and E.B have fixed values. Particularly, if E.B = 0, then we can always find a reference 
frame in which E = or B = (according as E^ — B"^ < ot > 0), ie. the field is purely 
electric or purely magnetic. In the following discussion we'll first find the transformation 
laws for the graphene charge carriers and then we analyse the consequences of application 
of Lorentz boost |116[ I117[ 1119] . For this, it will be convenient to write graphene equation 
in Lorentz - covariant form |119] : 

Y{p^-a,)^ = (71) 
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where the 2x27^ matrices satisfy the anticommutation relation [7'^,7'^]+ = g'^'^', g'^" is 
the metric -tensor. We can exphtly write the gamma-matrices as: 7*^ = cr^, 7^ = —iay, 
7^ = —i<Jx, where <Jx,<Jy, <Jz are Pauh matrices. Here we use the space-time notation for 
coordinates as = {vFt,xi,X2), momenta = h{ivp^dt, —idxi, —^9x2), and external field 
= («o, «i, 0.2) = (^~^^ 

The above EqjTlJis invariant under Lorentz transformation. The transformation laws can 
be obtained on the same lines as for the general (3 + 1) Dirac equation. The space-time 
coordinates, momenta and fields transform as: 



x"^ 



(72) 



where S{A) = exp d'^fiuil'^,!'^]) for A = exp^u). 

We consider the magnetic field modulations where the field strengths vary in x direction 
and are constant in y-direction and scalar potential V = V(x), then ai = in EqjTTj Writing 
the solutions as \E'(x, t) = ^'^^vVQ-'^Et/h^ obtain 



Ei){x) 



-ihvpaxdx + vpfTy {hky H — A{x)^ + V{x) il){x) (73) 

Now we apply two set of transformations [TTT] (A) a lorentz boost in y direction with 
rapidity z7r/2, followed by (B) a reflection about the y-axis. These two transformations 
together make a complex Lorentz boost as det{AA)det{AB) = 1. With this complex Lorentz 
boost applied the energy and momentum transforms as : 

E = —ihvpk'y 
ky = lE'/hvp 



And transformed Eq. 73 becomes 



E'ij\x) 



-ihvFCTxdx + o-y {f^vpk'y — iV {x)) H -iA{x) ip'{x) 



(74) 



On comparing 73 and 74 it can be clearly seen that it is the imaginary scalar potential 
which is now coupling with the y-component of momentum and the role of scalar potential 
is now played by the external vector potential (again imaginary). 

This section concludes our reviewing various aspects of ballistic electron transport in 
monolayer graphene through magnetic barriers, tuned by gate voltage and the associated 
optical analogy. In the subsequent section we shall discuss the electron transport in bilayer 
graphene in presence of such magnetic barriers. 
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VII. MAGNETIC BARRIERS IN BILAYER GRAPHENE 



A. Transmission through magnetic barriers 

In this section, we consider the effect of such inhomogenous magnetic fields for the case of 
unbiased bilayer graphene (BLG). Bilayer graphene is modelled as two coupled hexagonal 
lattices including inequivalent sites a, b and c, d in the bottom and top layers respectively. 
The present work discusses bilayer graphene in bernal stacking. In such an arrangement, the 
bottom and top layers are shifted with respect to each other in such a manner that c sites are 
exactly above the a sites and b sites are exactly in the middle of the hexagons of the bottom 
layer. In addition to intralayer interactions, there are interlayer interactions present in the 
system. The band structure is derived after taking into account two possible ways of A ^ C 
hopping |120j : via the dimer state or due to weak but direct a — c coupling. There are other 
weaker tunnelling processes [26] which are neglected in this tight-binding approximation 
which we consider. The tight-binding Hamiltonian |12U] derived in this way spans over the 
basis states that are located at a, b, c, d. The effective low energy Hamiltonian near the Fermi 
level |120] obtained within the k ■ p approximation from this tight binding Hamiltonian in 
this case is different from the corresponding one given in Eq.([5]) for monolayer graphene. 

The effective 2x2 Hamiltonian is given as 



The solutions of these equations are massive Dirac fermions since they have a quadratic 
dispersion and non-zero effective mass m. The applicability of this model Hamiltonian to 
study transport through a scalar potential was discussed in detail in |122j . In the presence 
of scalar potential V{x), under normal incidence ky = 0, pseudospin is still conserved, 
but the Hamiltonian is like that of a Schrodinger particle; i.e., H^fj = + V{x)I. 

The probability amplitude of corresponding charge carriers exponentially decays like usual 
Schrodinger particles inside the barrier leading to negligible transmission p ^[38lfT2T] . Draw- 
ing an analogy with optical transmission we call this phenomenon Klein reflection. 

Below we calculate the effect of magnetic barrier for the case of bilayer graphene (BLG). 
The electronic wavefunction depends on the spatial profile of the vector potential and the 
boundary conditions when an electron is incident from a region with i? = to a region where 
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B ^ 0. We show this by exphcitly calculating transmission by a transfer matrix approach 
[123] . jl24] . Chiral charge carriers in BLG obey: 



dt 



H 



Vi li et 

n+ 

et V2 n+ 

n ^2 



(75) 



Here, they are described by a 4-component spinor = (\[/^ v]>j^ v]>^ ^I/d)"^ a 4 x 4 

Hamiltonian in the presence of a magnetic barrier and an electrostatic potential. 11 = 
vf[Vx + i{py + 6v4/c)] with tij? = 10^ m/s and is the tunnel coupling between the two 
layers. Vi and V2 are the potentials at the two layers. The field profile taken here is 
B = BQ^x^ — (P)z. As pointed out earlier that such a sharp profile is routinely created 
in magnetic recording media where bit lengths can be varied over 10 — lOOnm range using 
suitable domain engineering [501 1125] . Since this length scale is much larger than the A-A 
lattice constant of 0.246nm in graphene but also much smaller than the Fermi wavelength 
\f = 50nm (at Ep = 17meV) of the incident electrons in bilayer graphene, the inhomogenous 
field pattern can be well approximated by the sharp edged profile proposed here. For bilayer 
graphene also, such sharp fields produced either by a patterned media layer or by a close- 
packed nanowire array can be easily integrated in the Si02 underlayer present in typical 
top- and bottom-gated graphene devices for Klein tunnelling studies The fields and 

voltages needed are within range of values already being used. 

We are interested in the stationary state solutions of the form \1/ = ijj{x,y)e~^^*^'^. Since 
[H,py] = 0, we can assume solutions like \I'a(x, y) = 0a(x)e*'^^^, for a = a,b,c,d. Substitut- 
ing this into Eq, 75, and using Ib = \/% and = ^ as units of length scale and energy 
scale respectively, to define the dimensionless quantities x — )■ ^ , 



V12 — , e'. and 



e — > we obtain the following set of coupled equations: 



i[d/dx - {ky + eAylB/hc)](j)b + t[(j)c = (e - 5)(j)a 
—i[d/dx + {ky + eAylB/hc)](f)a = (e' — S)4>b 

-i[d/dx + {ky + eAylB/hc)](j)d + e[(f)a = {e + S)(f)c 
—i[d/dx — {ky + eAylB/hc)](j)c = (e' + S)<Pd 



(76) 



57 



t 



o 



>> — 




-1 




+B 





+B 




B = 




B = 


B = 


Ila 


lib B = 


Region I 

Region II 


Region III 


Region I 




IRegion III 


(a) 


• 




(b) 


1 






100 
Y{mY) 



i: 



FIG. 25: Transmission, T Vs y and qy for various fields, (a) Barrier, (b) Barrier +well. A 
d = 50nm, Ep = 17meV. 



Here 



sgn{x)Bd/2 if|x| > ^ 



Bx \i\x\ < I 

fo with vq = {vi + V2)/2 = and 5 = {vi — f2)/2 



sgn{x) represents sign of x, e' = 

The solutions can be exphcitly obtained in the three regions a; < — | (region I), — 
a; < I (region II) and x > | (region III) as follows: 
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Region I(j = 1) and III(j = 3) are constant vector potential regions. We substitute 
(x) = 00 (0) 6*'^"=^' in Eq, 76 With this, the dispersion relation becomes: 



[-q{xY + (e' + 6Y][-q{xY + (e' - 5)^] = {e" - 6')t 



(78) 



Here, q{xY = ql + qy with qy{x) = ky + sgn{x)n-^. Also, ^ is the total magnetic flux 
through an area dls in units of the flux quantum $o = ~- The above Eq, 



7^ 



leads to 



propagating as well as evanescent wave solutions. The unbiased BLG (which is the case 

considered here), corresponds to 5 = 0, and vi^2 = f o = inregion I and III. Then the 
complete wave function in region I and III can be given as: 

^,{x,y)=Mo{x)e''^y{a, b, c, d,f (79) 



Mn(x) 



[qx - iqy{,x)]e^'i'''' -[qx + iqy{,x)]e ^''"'^ i[Kx - qy{x)]e ''"^ + gj;(x)]e''"'= 

-[qx + iqy{x)]e''^-'' [qx - iqy{x)]e-''^-'' i[Kx + qy{x)]e-''^'' -i[Kx - qy{x)]e''^'' 
For the wavefunction solutions in region II, where the magnetic field is finite, we introduce 
{x + ky) = z/y/2 and rewrite Eq, 76 as: 



iV2[d/dz - z/2](pb + e>c = (e' 
-iV2[d/dz + z/2](j)a = {e 



5)06 



-iV2[d/dz + z/2](Pa + e[(P, = (e' + 5)0, 
-iV2[d/dz - z/2](pc = (e' + 5)(pd 

With this, we obtain two coupled equations in (pa and 0c which looks like 

((e' + 6Y + 2{d^/dz^ - z^/A - 1/2)) 0, = (e' + 5)e>, 
((e' - 5f + 2{d^/dz^ - z^/A + 1/2)) 0, = (e' - 5)40, 



(80) 



(81) 



The above two equations can be solved to obtain a fourth order differential equation 
which can be decomposed in two second order differential equations. For (pa, we obtain : 

{d^/dz^ - z^/A + -f+/2){d^/dz^ - z^/A + 7_/2)0, = 0, 

7± = + 5' ± [(1 - 25e')' + (e'' - 5^^^'^ (82) 
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The above equation admits the solutions which are of the form of parabolic cylindrical 
functions [TO]. For the case of unbiased BLG the complete wavefunction solution in region 
II can be given in terms of a compact matrix as: 



Mb{x) 



elp-^Dp+^iiz) elp~Dp-_i{z) eip+Dp+_i{-z) Sip" Dp-_i{-z) 

e^Dp+{-z) e^Dp-^-z) 



stDr 



-2 V^J 



eie2Dp+^i{z) £i£2 -Dp-+i(2;) ele2Dp+^i{—z) er^eg-^p^+il 



Here D corresponds to parabolic cylindrical function of argument z and order p^ = 
(7± - l)/2. Also £i = ^ and = - and 7± = e'^ ± yiTe^. The current 

(Tx 



density expression is obtained as = vpi^^ 







ip. The transfer matrix through any 



combination of a scalar and vector potential can now be written in terms of transfer matrices 
Aisix) and AiQ{x) for regions with finite and zero magnetic field respectively. This can 
then be used to find the transmission. The ratio of current density in region III and the 
incident current density in region I gives the transmission probability as a function of the 
angle of incidence 6' = tan~^ — . 

^ qx 

In region I we parametrize q = {qx, Qy) in terms of incident energy, (e > 0) and the angle 

In region III similar 



of incidence 0' as {qx, ky — 



cos ( 



'eel + e"' sm i 



parametrization can be done in terms of the angle of emergence 9' to write {Qx, ky + ^) = 
{\/ee^ + e^cos^', A/ee'^ + e^ sin^'). The conservation of energy and constancy of the normal 
component of momentum then ensure 



sin 9' 



d 



+ sin( 



Ut[ + e2 

This equation explains that beyond a certain critical angle the electron beam and its 
cross-section (allowed values) is reduced with increasing barrier strength. This feature, 
along with chiral symmetry breaking, dictates the nature of transmission as is subsequently 
explained. 



In figure 25 is plotted transmittance T for a magnetic barrier ((a), left column) and for 
a barrier + well ((b), right column). In figure 1(a), the uppermost plot gives T for B = OT, 
clearly showing the region of perfect reflection around normal incidence symmetrically placed 
between two wings of resonant FabryPerot fringes. This is a generic feature of transmission 
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V{mV) V{mV) — 

FIG. 26: G for different magnetic barrier strengths, (a) Barrier, (b) Barrier +well. A = 0, d = 
50nm, Ep = 17meV. 

through such barriers |41j|123] . As the barrier strength increases to B = 0.2T, a transmission 
region develops between these two wings due to the effect of the magnetic barrier and 
the resulting transmission also becomes asymmetric. For higher fields, the disappearance 
of perfect Klein reflection at normal incidence is seen even more clearly. We may then 
conclude that in BLG, though the passage from electron to hole states at normal incidence 
is usually forbidden due to pseudospin conservation, even a weak magnetic field lifts this 
restriction. For a weak magnetic field, the coupling between electron and hole states is also 
weak; however, for a case of resonant scattering one can expect total transmission |126] The 
total angular range of transmission however shrinks in the presence of the magnetic barrier 
since, beyond a certain incident angle, all electron waves suffer total internal reflection. In 
figure 1(b) is plotted T through a barrier + well, consisting of two magnetic barriers equal 
in magnitude and opposite in direction such that the total flux through the region vanishes. 
Now, the incident and the transmitted wavevectors are parallel to each other. As a result, 
the net rotation of the pseudospinor due to the inhomogeneous field vanishes and Klein 
reflection at normal incidence is restored. The FabryPerot fringes bend due to asymmetric 
transmission at a given V. 

B. Conductance 

We now study the effect of the above transmission on conductance at very low temper- 
atures and for energies close to the Fermi energy in the linear transport regime. To a good 
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approximation, the dimensionless conductance can be written as [73] 

TT 

G= rf0r(E,0)cos0 (83) 



Fig, 26 plots conductance G as a function voltage. This shows that the gaps between the 
conductance maxima and minima get reduced in the presence of a magnetic barrier. The 
barrier reflects electrons incident upon it beyond a critical angle and the angular range of 
transmission shrinks with increasing barrier strength. As a result, the absolute value of con- 
ductance maxima comes down. Also the asymmetry, namely the way G changes to the right 
and left side of such maxima, reduces in the presence of magnetic barrier. The behaviour 
of the conductance is strongly dependent on transmission resonances. To explain this be- 
haviour, therefore we shall now model these numerically computed transmission resonances 
using a the Breit Wigner model |127] 

Such a form is valid near the transmission resonance. The above model has been widely 
used to model the transport due to resonant tunnelling in semiconductor structure [T]. 
Recently this model has been used to understand the cloaking phenomena in bilayer graphene 
|128] . Here e'^{ky),V{ky) are respectively the momentum dependent resonant energy and the 
resonance width. 



Following the expression (78) in the absence of a magnetic field, the resonance energy 
^*{ky) and the resonance width T{ky) can be assumed to have a quadratic dependence on the 
wavevector ky |128j . Thus e'^{ky) ^ e'^(0) + a/c^ and T{ky) = f3ky. The quadratic dependence 
of the resonant energy and the transmission width explains why the Fabry-Perot fringes of 
transmission resonance depicted in the left column of Fig. [27] are of convex shape. The 
asymmetry in the form of this fringes reflects themselves in the plot of the conductance 
again depicted in the left lower corner of Fig. [27} as one goes from the conductance peak 
to the left and right along the e' axis (effectively the voltage axis). In order to analyse the 
conductance peak, we evaluate the quantity 6G{e') = G(e') — G^e'^^), Here G^e'^^) repre- 
sents the conductance value at the peak, and G(e') represents the conductance value in its 
neighbourhood. Then, 



SG{e') = Go J {\t{e',ky)\'-\t{e:„ky)\'^ dky 
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FIG. 27: Explanation for the conductance: The upper row consists of zoomed in portion of the 



plot given in Fig. 25 and the lower row consists of the zoomed in portion of Fig. 26 



On substituting the quadratic dependencies of reseonance energy and resonance width in 
the integral, 



6G{e') = / dk 



which upon solving gives: 



cos[^ + ^sgn(Se')] \Se'\ , , , 

^g(^0 = -govr ^(\^l,)3/2 « = «//5 (85) 

Clearly ^^^^ shows the sqaure root sigularities at the conductance peak and the sgn{6e') 
term appearing in the expression introduces the asymmetry on the two sides of resonance 
energy. 

In the presence of a finite magnetic barrier, the small ky behaviour of such resonances 
can be obtained by using a semiclassical argument. We replace ky by its covariant form, 
namely ky + 2tt^ . This yields e'^{ky) ^ e'^(O) + Ad'^ where A is a constant. With increasing 

e 

strength of the magnetic barrier, e'^(O) is shifted by the amount of the flux inserted in 
the barrier regime and this shift is independent of ky. The resonance width also follows a 
similar behaviour. The right column of figure 3 depicts this situation. The field-induced 
resonance leads to a central fringe symmetric about the voltage axis, unlike the upper and 
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FIG. 28: Selective blocking of modes in trilayer graphene. E and k are in arbitrary units. The 
scalar potential V blocks massive modes while the massless mode is passed through. The massless 
mode is weakly transmitted through the magnetic barrier B. 

lower branches, which are asymmetric. This is again seen in the conductance plot as a right- 
left asymmetry around the conductance peak, which is now reduced. The reduction of the 
conductance peak can be attributed to the shrinking of the angular range of transmission due 
to the magnetic barrier, whereas the rounding of the conductance peak and the consequent 
reduction in the left and right asymmetry around this peak is an indicator of the reversal 
of Klein reflection. Since conductance is a physically measurable quantity, this provides a 
direct way to experimentally test our results. 

C. Generalisation beyond bilayer 

The above findings also extend to multilayer graphene. For N- layer graphene with Bernal 
{ABAB ■ ■ ■) stacking The dispersion relation is given by |129] 

where r = 1, 2, 3, ■ ■ ■ corresponds to the N modes of the energy spectrum, vp is the Fermi 
velocity. For odd N, r = is massless and it's energy is given by 

Et^ =±vp\p\ (86) 

2 't^ 
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For all other — 1 cases the dispersion relation takes the form E"^ + 2£'t cos(-^^^) = ti^lpp. 



In the limit \E\ « the low energy spectrum for the r-th massive mode is given as 

2 

p rvr 
Erv = , if Qcosf ) < 

rvr 

,if etcosf— -) < 



2m/ + 

where we define rrir = ^"^^ ^°'^(iv+i)l ^j-^ effective mass of the r-th massive mode. Such band 



structure has been plotted in fig. 28 for trilayer graphene. Now the transport electrons near 
the Fermi level will be a mixture of eigenmodes, one of which is a massless Dirac mode 
and the others are massive modes with different effective masses and also have a different 
chirality. These arguments extend to > 3 also. An extension of the semiclassical argument 
presented earlier therefore suggests the following possibility in the case of such multilayer 
graphene. When such a system is exposed to a scalar potential barrier, a normally incident 
massless mode will undergo Klein tunnelling, whereas in the case of a massive mode there 
will be exponentially decaying transmission. In contrast, a magnetic barrier will be highly 
reflective for a massless mode but will transmit normally incident electrons for massive 
modes. A typical situation has been depicted in fig|28| Since the low energy band structure 
decomposes into one massless and A^ — 1 massive modes for N odd and only N massive modes 
for N even [129J, our results qualitatively suggest that a magnetic barrier can be used to 
selectively allow or partially filter out a mode. 

In Fig j29] we compare the modified bandstructure obtained with a periodic arrangement 
of DMVP in monolayer and bilayer graphene. One dimensional (ID) periodic modulation 
of the chemical potential or the electric field perpendicular to the layers was addressed in 
Ref. |130] . For a periodic stratified media of several DMVP barriers, forbidden region 
appear in the bandstructure where the propagation of Bloch waves is not allowed. Th size 
of these forbidden regions increase with increasing strength of magnetic field. A comparison 



between the plots in Fig, 29 (a),(c) and (b),(d) shows very clearly the conical band structure 
of monolayer graphene and parabolic bands in bilayer graphene. 

This issue of anti-Klein tunnelling effect in bilayer graphene p-n-p junction and the effect 
of magnetic field on it, has also been addressed Ref. |128j . Here they describe it in terms 
of the bound states being completely decoupled from outside continuum states, thus calling 
it an electronic cloaking effect, and the barrier acting as "cloak". The effect is complete at 
normal incidence, leading to "Klein reflection" or "anti-Klein tunneling". This is explained 
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FIG. 29: The band structure modification of monolayer and bilayer graphene in presence of mag- 
netic barrier, (a) and (c) plots the bandstructure of monolayer graphene in presence of periodic 
arrangement of magnetic barriers with two different strengths where as (b) and (d) plots the 
corresponding bandstructure of bilayer graphene. 

as follows: At normal incidence the effective equation that describes the motion of charge 



carriers in graphene , (as can be seen from Eq 75), can be analysed in the eigen basis of a^'- 

±[V-E]]i,^{x) = Q (87) 



2m dx"^ 



Clearly, the solutions for ipj^ corresponds to the solutions outside the barrier. And the 
solution for (for which the equation appears with a negative sign) are bound state 
solutions, and corresponds to the solutions inside the barrier. Thus these two kinds of 
solutions are completely decoupled leading to zero transmission despite of the states available 
inside the barrier. 



In section |VII[ our discussion on Klein reflection in bilayer graphene was based on the 
analysis of Fabry-Perot like transmission resonances developed from inhomogenous magnetic 
barriers in the presence of voltage. Here the oscillatory transmission due to Fabry-Perot like 
arises due to interference of forward and backward moving propagating waves which undergo 
multiple reflection between two interfaces. In addition to these FP like resonances a new 
kind of transmission resonances have also been seen for the case of gapped bilayer graphene 
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[131] ■ In their recent work, R. Nandkishore et. al. proposed an entirely different approach to 
reahze oscillatory transmission which involves only forward propagating waves and a single 
interface and thus is different from the FP resonances. They show that Zener transmission 
in a p-n junction in gapped bilayer graphene exhibits common path interference taking 
place under the tunnel barrier, leading to transmission that oscillates as a function of gate- 
tunable bandgap.The origin of this oscillatory behaviour lies in the fact that in the presence 
of p-n junction, bilayer graphene employs the propagating as well as evanescent wavevector 
solutions. The simultaneous existence of these two- leads to common path interference which 
gives rise to oscillatory transmission. These oscillations in transmission manifest themselves 
through negative differential conductivity in the I — V characteristic. Interesting works were 
also done on pseudo spinotronics in biased bilayer graphene |132[ llSSj which we shall not 
discuss in this review. Other recent developments made on these lines include the work by 
Campos et. al. |134] . where they have realized giant conductance oscillations in ballistic 
trilayer graphene Fabry-Perot interferometers, which result from phase coherent transport 
through resonant bound states beneath an electrostatic barrier. 

VIII. CONCLUSION 

A large body of work has followed since the proposal of using inhomogeneous magnetic 
fields in graphene and understanding the transport in an analogy with light beam propa- 
gation in optics- termed as electron optics. The similarities between polarization states of 
light and ballistic charge carriers in graphene has been investigated in detail in |135l 1137] . 
It is demonstrated theoretically |136] that regions with different electrostatic potentials in 
graphene can have the same transmission probability as the transmittance of layered struc- 
tures in optics illuminated with normally incident TE waves if the incidence angle of electrons 
in graphene and the width of the gated region are appropriately chosen. This quantitative 
analogy can be useful for designing optical structures that correspond to graphene-based 
devices. Then there are proposals of electron waveguides created by electrostatic potential 
or by real magnetic barriers [S^ I138j . Here the control the of the quasiparticle flux in a 
graphene-based waveguide was theoretically investigated. Guided modes occuring in the 
negative-zero-positive index metamaterial waveguide [141 j and in symmetric velocity bar- 
rier [140 ], have also been investigated. It is shown that magnetically induced waveguide 
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in graphene leads to strong confinement of Dirac fermions, regardless of its edge termina- 
tions |142] .The properties of unidirectional snake states have also been investigated in such 
waveguides |139] . It was shown that for certain magnetic field profile, two spatially separated 
counterpropagating snake states are formed, leading to conductance quantization which is 
insensitive to backscattering by impurities or irregularities of the magnetic field. 

It may be mentioned there are many other interesting articles in the relevant field and 
because of space restriction it is almost impossible to refer to all of them which we regret. We 
would also like to direct the readers some excellent reviews on electron transport in graphene 
where some of the related issues was discussed in somewhat different context. In ref. |144j 
the transport properties of graphene that includes theoretical and experimental work was 
reviewed. In ref. |145] the charge and spin transport in mesoscopic graphene structures was 
reviewd in detail. In ref. |146j the density and temperature-dependent charge transport in 
doped or gated graphene devices was reviewed. In ref. [23l [23] the Klein tunnelling through 
single and multiple potential barriers was reviewed. Properties of graphene in presence of 
strong homegenous magnetic field was reviewed in ref. |147j . Electronic properties of bilayer 
graphene was reviewed in detail in ref. |148] . 

To summarize. We have reviewed several aspects of charge carrier transport in graphene 
through inhomogenous magnetic fields viz. magnetic barriers in detail. We reviewed ex- 
perimental progress in designing such magnetic barriers at various length scales. Such bar- 
riers suppress Klein tunnelling, thereby achieving confinement in graphene, which can be 
seen through strong supression of transmission of electrons. We have shown that transport 
through MVP barriers can be understood in terms of propagation of light through periodic 
stratified media. This analogy can partially be attributed to the fact that equation describ- 
ing Dirac like charge carriers in graphene and the Maxwell equations |143] are both linear 
wave equations. We show that such transport properties through a singular magnetic barrier 
can be much better controlled by the additional application of an electrostatic voltage. A 
detailed optical analogy for transport through single and multiple EMVP barriers arranged 
periodically was obtained highlighting optical analogues of phenomena such as TIR for pos- 
itive and negative refraction and a Quantum Goos-Hanchen shift. One significant feature of 
such transport in this regime of the simultaneous application of highly inhomogenous mag- 
netic fields along with electrostatic potentials, is highly anisotropic and strongly dependent 
on the sign of the voltage indicating possible device applications. We discuss coUimation, 
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emergence of extra dirac points in this context. We also reviewed transport through such 
magnetic barrier through unbiased graphene bilayer, describe reversal of Klein tunnelling, 
anticloaking, common path interference etc. in such systems. The wide range of optical 
analogues to charge transport in graphene may lead to practical devices such as Bragg re- 
flectors, resonantors, waveguides etc. We hope this detailed review will augment theoretical 
and experimental research in this field. 
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